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ABSTRACT 


An experimental investigation is conducted to examine the simultaneous 
motion of three phases in porous media. Limited research has been done in this 
area; therefore, the work is exploratory in nature. Experimental apparatus is con- 
structed to enable the flow of three immiscible fluids through a one-dimensional, 
horizontal glass bead packed core. The fluids used are water, benzyl alcohol and 
decane, representing water, oil, and gas phases respectively. The apparatus allows 
the core to be brought to any combination of saturations of the three phases through 


steady-state and unsteady-state (dynamic displacement) experiments. 


The main objective of the research is to map the relative permeabilities of the 
fluids in the interior region of the ternary diagram through a series of dynamic dis- 
placement water floods. Fluid relative permeabilities and saturations during the 
dynamic displacement experiments are calculated from the three-phase extension of 
the two-phase dynamic displacement theory of Welge’s permeability ratios, and 
Johnson, Bossler and Naumann’s two-phase method of determining individual rela- 
tive permeabilities. Steady-state relative permeabilities are determined from 


Darcy’s law. 


The interior region of the ternary diagram cannot be mapped by a dynamic 
water flood originating along the residual water saturation line. The piston-like dis- 
placement of decane and the development of a benzyl alcohol (oil) bank create a 
Buckley-Leverett saturation jump which follows the residual water saturation line 
down to the dynamic residual decane saturation line. The saturation trajectory then 
follows the dynamic residual decane saturation line down to residual benzyl] alcohol 


and decane saturations in the presence of water, bypassing the central region of the 





ternary diagram. A difference in residual decane saturations is seen between the 


steady-state and the dynamic displacement floods. 


Exact graphical matching of the core-end saturation change is seen between 
the displaced decane and the oil bank. This observation lends credence to the 


extension of the two-phase dynamic displacement theory to three phases. 


Further investigation will be required to examine the mechanism behind the 
development of the oil bank and the difference in the residual decane saturations of 


the two types of floods. 
Implications of the discovered oil bank phenomena include examining alterna- 
tive methods for mapping the interior of the ternary diagram and possible applica- 


tions to oil recovery. 
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Chapter 1 


INTRODUCTION 


The objective of this research is to further our understanding of the simultane- 
ous motion of three phases in porous media. This is to include hysteretical effects 
and to specifically quantify and qualify immiscible three-phase relative permeabili- 


ties. 


Multi-phase flow in porous media occurs in many engineering fields, such as 
ground water hydrology, hazardous waste disposal sites, surface spills, and 
petroleum hydrocarbon reservoirs. As the upper crust of the earth is being taxed 
more and more with sophisticated processes, a need for a better understanding of 
transport phenomena in porous media is required. We have a responsibility to 
preserve the functionality of fluid bearing formations, and to be able to accurately 
predict and optimize the outcome of prospective production or injection projects. In 
the case of contamination, we must have the ability to design a cleanup process, 
minimizing damage and danger. 

One of the key parameters used to describe the motion of various fluids in 
porous media is relative permeability. The relative permeability parameter is used 
to determine the energy needed to move the fluids through the porous media, and 
the relative flow rates of each phase present. Knowledge of the permeability of a 
formation rock to a given fluid is an essential input for numerical models used to 
predict and optimize aquifer and reservoir performance, and for the designing of 


efficient and effective mobilization or utilization systems. 





Several examples of three-phase flow in subsurface formations are: 


е А Surface spill of a non-aqueous solution, not miscible with water, that seeps 
into the ground under gravity. This non-aqueous solution then mixes with rain 
water and air, to form a three-phase mixture flowing simultaneously in a 


porous medium. 


e A non-aqueous fluid might be introduced into a disposal well that penetrates 
an aquifer containing air or a non-condensible gas such as carbon dioxide. 


Again, the result is a three-phase flow. 


e In hydrocarbon production, simultaneous flow of oil, gas and water is com- 
mon, under both primary or secondary recovery. Three-phase flow is also 


encountered in enhanced recovery processes. 


e Three phases can be created in the subsurface in an effort to control the move- 
ment of a contaminant toward a set of producing wells, through the use of 


some displacement process. 


The relative permeability to a given phase is a function of that phase’s satura- 
tion and possibly, a function of other phase saturations as well. The effective per- 
meability is one of the major non-linear factors in the mass transport equations that 
are used in both analytical and numerical reservoir models. Due to the highly non- 
linear nature of the multiphase flow problem, relative permeabilities must be deter- 
mined experimentally. The results and trends from this exploratory experimental 
research have application to multiple disciplines. Therefore, in this work a very 
idealized system consisting of a one dimensional core, simple fluids, and relatively 


homogeneous porous media are used. 





Chapter 2 


REVIEW OF THE THEORY AND LITERATURE 


One-dimensional single-phase flow through porous media is described by 
Darcy’s equation (Darcy 1856), shown as: 


паче (2.1) 


Darcy’s law applied to all phases (capillary pressure neglected) is: 


END 
ве, (2.2) 


Vi 


where k; is the effective permeability of phase 1. 


Relative permeability is defined as: 


k; 


ki = (2.3) 


КЪ 





where k,; is the relative permeability of phase 1. 


2.4 Steady-State and Dynamic Relative Permeability Measuring Techniques 


There are two major methods involving continuous injection of fluids for 
determining relative permeabilities: the steady-state method and the unsteady-state 
method. In the Penn State steady-state injection technique, a fixed fractional flow 
of fluids is injected into the core, until the same fractional flow of fluids is pro- 
duced from the core. At this point it is assumed that steady-state is achieved. The 
pressure drop across the core is measured and Equation (2.2) is used to determine 
the effective permeabilities followed by Equation (2.3) to determine the relative per- 


meabilities of each of the injected phases. 





The unsteady-state method of relative permeability determination (dynamic dis- 
placement) involves the injection of (usually) one fluid phase into an iso-saturated 
core. The fluids recovered from the core, and the pressure drop along the core, are 
measured as functions of the pore volumes injected. From this data the core-end 
saturations and the associated fluid relative permeabilities are determined. The 
specific method of determining these values is derived by extending, from two- 
phases to three-phases, the two-phase dynamic displacement theory of Welge’s per- 
meability ratios with Johnson, Bossler and Naumann’s (JBN) two-phase method of 
determining individual relative permeabilities (Welge 1950; Johnson et al. 1959). 
The actual derivation and extension of the theory is presented by Grader and 


O’Meara (1988) and is summarized in section 2.3. 


2.2 History of Three-Phase Relative Permeability Determination 


A large number of papers relating to the determination of relative permeabili- 
ties in porous media is found in the literature. Honarpour et al. (1986) presented a 


review of the subject. 


Leverett and Lewis (1941) presented a study of three-phase relative permeabil- 
ities using the steady-state method. Their conclusion is that the water relative per- 
meability is solely a function of the water saturation in an unconsolidated sand 
pack. 

The steady-state method for determining three-phase relative permeabilities 
was used by Caudle et al. (1951), Corey et al. (1956), Snell (1962), and Oak et al. 
(1988). Caudle et al. (1951) found that for a consolidated core, each relative per- 


meability depends on the saturation of the three phases. 


The few papers listed above used the steady-state method for determining 





three-phase relative permeabilities. Measurements of relative permeabilities using 
the unsteady-state method are presented by Sarem (1966), Donaldson and Dean 
(1966), and Saraf et al. (1983). Sarem (1966) used a Welge-like method which 
assumes that each phase relative permeability depends only on its own saturation. 
As shown by Grader and O’Meara (1988) this assumption is not needed in the 
development of the extension of the Welge method from two phases to three 


phases. 


Donaldson and Dean (1966) used a variation of Sarem’s method (Sarem 1966) 
in that the relative permeabilities depend on the saturation of the three phases in the 
core. For the three-phase case, they used the Welge method (Welge 1950) by 
replacing the fractional flow of oil, by the sum of the fractional flow of oil and 
water. They do not, however, present a derivation or the assumptions made for the 
application of the Welge method to three phases. In the paper the three-phase rela- 
tive permeabilities for Berea sandstone and Arbuckle limestone are presented. 
These two rock types show high residual water and oil saturations, limiting the size 
of the three-phase region on the ternary diagram available for mapping. Saraf et al. 


(1983) also make use of an extension of the Welge method. 


Summary tables with some of the three-phase relative permeability work 
described in the literature are presented by Honarpour et al. (1986), Oak et al. 
(1988), and Donaldson and Kayser (1981). 


The difficulty in determining three-phase relative permeabilities using either 
the steady-state or unsteady-state methods justifies the need for a three-phase rela- 
tive permeability model which uses only two-phase relative permeability data. 
Baker (1988) presents a comprehensive summary of the various three-phase relative 


permeability models. 





2.3 Extension of the Welge/JBN Theory to Three Phases 


The Welge/JBN method discussed in this section is based on the work of 
Grader and O’Meara (1988) which extends the Welge/JBN theory from two phases 
to three phases. The theory is simplified to the non-gravity case. In the following 
derivation, the standard simplifying assumptions are made: incompressible, one- 
dimensional, isothermal, and immiscible three-phase flow with no capillary pressure 
effects. The material balance equation and Darcy’s Law equation describing the 


velocity of fluid flow in a porous medium are: 


05; Оу; 


ф E E ЕЕ = 0 (2.4) 
д 
y=-4 3 (2.5) 


where j=1,2,3, S¡+5S2+S3=1, A¡=k;,/H;, and the rest of the terms are 
defined in the Nomenclature. By summing the three material balance equations we 
can show that v, the total fluid velocity, is independent of x, and is therefore, at 
most only a function of time (we will show later in the Chapter 7 how the variable 
rate does not affect the validity of the proposed solution). Summation of Darcy’s 
Law equations for the three fluids provides an expression for the pressure gradient 


in terms of v, which can then be substituted into Eq. (2.2) to obtain: 


У; = У D (5152) (2.6) 
where 
^j 
А, = А + А + А (2.8) 


Notice that, for a given choice of fluids, f; is a function only of saturation. Also, 


note that there are only two independent volume equations for two phases, as the 








third phase volumetrics are complimentary, and dependent. Substituting Equation 
(2.6) into Equation (2.4) yields the following differential equation (neglecting grav- 
ity effects): 


ES E Ox 


For suitable initial and boundary conditions (iso-saturation at initial time, and a 


às, f әб | аз, ы” X 
MES tivos, = or |- |І, (2.9) 


constant fractional flow at the injection face), the solution to Equation (2.9) is self 


Q 


similar. A natural choice for the similarity variable, n, is = where X denotes the 


dimensionless position scaled by 1, the length of the core and Q denotes the cumu- 


lative pore volume of total fluid injected: 


t 
‘joe [ «dr (2.10) 
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Substitution of ņ into Eq. (2.9) yields: 
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Letting S? denote the initial saturation of phase j in the porous medium and 5) 
denote its initial saturation at the inlet end of the core, we obtain the initial condi- 
tion: 

SUIS at 7 = 0 (2.13) 
and the boundary condition: 
5, = 5) а Э е (2.14) 

Тһе Яша mobilities and consequently, the permeabilities are calculated from 


the following equation that is obtained from Equation (2.5): 





р (2.15) 


where 





Y == (2.16) 
дх 

Ф=- (2.17) 
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V; and @ are functions of only the fluid saturations and, thereby, n. Our task then 
reduces to experimentally determining the functional dependency of S;, yj, and P 
on Tl. We assume that the conditions for self-similarity prevail. Then from Equa- 
tions (2.6) and (2.16) we obtain: 
W; = fi (S1, S2) (2.18) 
Substitution of Equation (2.6) into Equation (2.5) shows that ® also depends only 
on fluid saturations. Therefore, for a self-similar system, Sj, y; and ® each depend 
solely on 1. 
Let L; represent the cumulative pore volume of fluid j collected at the outflow 
end (X = 1) of the porous medium. Then y; is simply the derivative of the L; vs. 


Q curve where n = Q, the fractional flow of fluid j: 


dL; 

— а. 
= ое fon a= © (2.19) 
' dQ dQ 

dt 


The average saturation of fluid j between the injection point (x = 0) and x is: 


с = (2.20) 


У X 
Upon differentiating Equation (2.20) we obtain: 














See 021) 
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Equation (2.21) relates the saturation at a particular point to the average saturation 
between the point and the injection end of the core. Since S; and, therefore, S; 
depends only on 1 this equation can be transformed in terms of the pore volumes 
injected as: 
Е 15 
Sj = Sj- Fm (2.22) 


The average saturations in the core, S; can be obtained via mass balance: 
S; = SP + y; Q = L; (2.23) 


where yi is the fraction of fluid j in the inflow stream. On substituting Equation 


(2.23) into Equation (2.21) for n = Q we obtain: 


е уй 2.24 
E uM (2.24) 
ог: 


Equations (2.24) and (2.25) describe the Welge method for determining Ше 
saturation at the down stream end of the core. Johnson et al. (1959) showed that the 
pressure gradient at the end of the core is needed in order to determine the indivi- 
dual relative permeabilities at the saturation values. Since @ is also a function of 


n, it too can be obtained from the measurement of its average: 


= Q4 d® 
Ф-Ф-1-- 2.26 
ШЕГЕ (2.26) 


At the end of the core X = 1, n = Q, the pore volumes injected. Hence, Equation 


(2.26) simplifies for pressure drops measured between the two ends of the core: 





2-2-055 (2.27) 


Finally, Equation (2.15) is simplified for the zero gravity case and for using the 


Saturation and pressure gradients at the downstream end of the core: 


J 
Where f; is the fractional flow of phase j, and 9 is the normalized pressure gradient 
at the end of the core. The application of the JBN method for determining the indi- 
vidual relative permeabilities can be used for a three-phase system, in the same way 


it is used for a two-phase system. 


2.4 Analysis Procedure 


The analysis procedure for determining the individual phase relative permeabil- 
ities and their corresponding saturations is similar to the two-phase procedure. This 
procedure was presented for three phases by Grader and O'Meara (1988) and is 
summarized here. The recovery of the three phases and the overall pressure drop 
between the injection and the production ends of the core must be measured as a 
function of time. The saturations and permeabilities are computed for the condi- 
tions at the end of the core. Figure 2.1 presents two synthetic recovery curves for 
two of the three phases in the system. The solid dots represent breakthrough. The 
Welge tangent construction method is described in this figure, similar to the method 
presented by Jones and Roszelle (1978). For a certain value of pore volume of fluid 
injected, the recovery of a given phase at the intersection of the tangent and the 
recovery axis is the change in the saturation of that phase at the end of the core. 
Note that there is no mathematical reason preventing this saturation change from 
being negative, indicating an increase in the saturation of the phase. The slope of 


the recovery curve with respect to the pore volume injected is the fractional flow 
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Figure 2.1: | Recovery of Two Phases Displaced by a Third Phase as a Function 
of PVI. 





value for a phase. We use the tangent construction method for two phases, using the 
individual recovery curves. The saturation of the third phase is determined using the 
material balance equation as well as its fractional flow value. Hence, Equation 2.24 


is used to evaluate 5; at the end of the core. 


Figure 2.2 presents the normalized pressure drop as a function of pore volume 
injected (the pressure drop is divided by the flow rate). In some cases discussed in 
the Chapter 7, the instantaneous flow rate was used, as the overall injection rate was 
not constant. For a given value of pore volume injected, the tangent construction 
method proposed by Johnson et al. (1959) (the JBN method) is applied. The inter- 
section point of the tangent with the normalized pressure drop axis yields the value 
of the spatial gradient of the normalized pressure drop at the downstream end of the 
core. This local gradient at the downstream end of the core is then used in con- 
junction with Equation (2.28) and the slope of the recovery curve to determine the 


relative permeabilities of the three phases: 


f lcore H ) 





АЕ (2.29) 
Kabs Acore 9 | 
Note that by using the definition of the similarity variable, 
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This is clearly shown in Figure 2.2. 
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Figure 2.2: Normalized Pressure Drop as a Function of PVI. 





Chapter 3 


PROBLEM STATEMENT 


The purpose of this research is to investigate the simultaneous motion of three 
phases in porous media. A series of laboratory experiments are conducted in a 
highly idealized environment, using; three immiscible fluids of similar densities and 
interfacial tensions, and a horizontal one-dimensional glass bead pack. The objec- 
tive of this work is to quantify and qualify immiscible three-phase relative per- 
meabilities. This is achieved by performing experiments to determine three-phase 
relative permeabilities using steady-state and non-steady-state methods, and by 
applying the frontal advance three-phase theory in the mapping of the ternary 
three-phase diagram. 
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Chapter 4 


EXPERIMENTAL APPARATUS 


4.1 Apparatus Capabilities and Overview 


The apparatus allows the performance of both steady-state and unsteady-state 
multi-phase flow experiments. The possible types of experiments include: 


a. one-phase injection 

b. two-phase Penn-State steady-state 

c. three-phase Penn-State steady-state 

d. two-phase Welge dynamic classical displacement 


e. three-phase Welge dynamic classical displacement 
Specific core saturations can be obtained through these types of experiments. Fig- 
ures 4.1 through 4.7 show the apparatus’ capabilities through the relation of the 
injected phase(s), initial core saturation, and resulting core saturation. The initials 
of the water, decane, benzyl alcohol used in the apparatus represent the 100% 


saturation points of the fluids on the ternary diagrams. 


Single-phase injections and saturations are shown by figure 4.1. The open dots 


are the injected phase and the resulting saturation of the core. 


Two-phase final saturations are possible through the injection scenarios shown 
in figures 4.2 through 4.4. Single-phase flooding into an existing single-phase 
saturation is shown in figure 4.2. The open dots are the existing phase saturations at 


the beginning of the flood and the solid dots represent the injected single-phase and 








Figure 4.1: 


Flooding Capabilities of the Apparatus; Single-Phase Injection 





Нриге 4.2: 





Flooding Capabilities of the Apparatus; One-Phase into One-Phase 





Figure 4.3: 
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Flooding Capabilities of the Apparatus; One-Phase into Two-Phases 





19 





Figure 4.4: Flooding Capabilities of the Apparatus; Two-Phases into Two-Phases 





Figure 4.5: Flooding Capabilities of the Apparatus; One-Phase into Three-Phases 
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Figure 4.6: Flooding Capabilities of the Apparatus; Two-Phases into Three- 
Phases 


Figure 4.7: Flooding Capabilities of the Apparatus; Three-Phases into Three- 
Phases 





the resulting two-phase saturations. Figure 4.3 diagrams the possibilities of a 
single-phase, shown by a closed dot, being injected into a two-phase saturated core 
and the resulting two-phase saturation, shown by an open dot. Two-phase satura- 
tions can be changed by injecting different fractional flows of the two fluids. Fig- 
ure 4.4 shows the two-phase saturation (open dot), being changed to a new two- 


phase saturation, (closed dot). 


Three-phase saturation paths are shown on figures 4.5 through 4.7. In figure 
4.5, a three-phase saturation, shown by an open dot can be achieved by injection of 
a single-phase, shown by a closed dot, into an existing two-phase or three-phase 
saturation. Three-phase saturations can be achieved or altered as shown in figure 
4.6 by the injection of a two-phase mixture into an existing two-phase or three- 
phase saturation. Three-phase saturations can be altered by the injection of a new 


three-phase mixture as shown in figure 4.7. 


A schematic representation of the experimenta! equipment used is given in 
Figure 4.8. Tubing is represented by thin lines. Valves are represented by dots. 
Important equipment is outlined by a bold line and described below. The fluid 
pumps; water (1), decane (2), and benzyl alcohol (3), draw fluid from either a cen- 
tral collection reservoir (4) or their individual fluid storage flasks. Each pump is 
flanked by a pressure surge suppressor to its left and its fluid storage flask to the 
right. The fluid injected into the system is kept below 20 psig by the use of a 
pressure relief valve (6). When pressure testing the system, the pressure is moni- 
tored up to 200 psig on the system pressure gauge (7). The pressure drop across 
the core (8), is measured with a digital pressure transmitter (9). Prior to its first 
flooding, the core is evacuated by a vacuum pump (10). Fluids leaving the core are 


returned to the system reservoir (4) or collected in a fraction collector (5). The 
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apparatus can be broken down into four major portions: 1. injection system, 2. 
porous medium/core holder, 3. collection system, and 4. pressure monitoring sys- 


tem. A detailed description of each portion is presented in this chapter. 


4.2 Fluids 


The three phases in this experimental work are represented by three immiscible 
liquids. The three liquids used are distilled water, benzyl alcohol (Phenyl 
Methanol), and decane. Distilled water represents a water phase, benzyl alcohol 
represents an oil phase, and decane represents a gas phase. The terms benzyl 
alcohol and oil are used interchangeably, as is decane and gas. All fluid storage 
vessels in the system contain all three liquids. Each pump has its own 2000 mL 
glass storage flask. Each pump storage flask contains approximately 1000 mL of its 
specified liquid and 500 mL of each of the other two liquids. The system reservoir 
consists of a 250 mL graduated cylinder and contains varying amounts of each 


liquid depending on the stage of the experiment currently being conducted. 


Рпог to each run the liquids in the storage vessels are circulated to allow 
presaturation with each other. Presaturation avoids interphase mass transfer during 


the run, and preserves the mass balance. 


The physical properties used for these three liquids are contained in Table 4.1. 


Table 4.1 
Summary of Fluid Properties 


p [g/ml] | p [g/ml] | p [g/ml] 
19°С 20° с 29956 


decane А 0.739 0.735 0.731 
water 1.124 1.002 1.042 1.001 
benzyl alcohol 4.943 4.217 0.992 1.026 1.022 
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The benzyl alcohol value for density at 15° c is not consistent with the tem- 
perature and is questionable. This has no effect as the experiments were conducted 


at about 23° c in the range where density values are correct. 


4.3 Injection System 


Three constant rate pumps are used to inject the three different phases, one 
pump for each phase. The pumps used are Beckman Model 110B Solvent Delivery 
Modules. These are fluid metering pumps designed primarily for liquid chroma- 
tography applications. Each module contains a single-piston pumping mechanism, 
drive motor, and electronic circuitry in one cabinet. Each pump can deliver from 
0.1 to 9.9 mL/min at pressures up to 6000 psi. A digital potentiometer allows 
selection of specified flow rates in increments of 0.1 mL/min. A comparison of 
actual versus specified flow rates is contained in the error analysis section, Chapter 
6, and each pump flow rate is calculated before each experiment. The pumps can be 
manually operated or remotely controlled by a suitable system controller. In this 


research project the pumps are operated in the manual mode. 


A valve system was constructed to allow the injection of one, two, or three 
phases. This valve system is mounted on a 1-foot by 4-foot piece of 1/4-inch plexi- 
glass sheeting to give the system rigidity. Each pump has a three-way inlet valve 
that allows the Ba to draw liquid from the pump’s storage flask or from the sys- 
tem reservoir (see figure 4.8). Downstream from each pump’s outlet there ıs a 
three-way pump outlet valve that allows the liquid to be returned to the pump’s 
storage flask or directed to the liquid selection valves leading to the core (see figure 
4.8) . These liquid selection valves determine if the liquid is to be mixed with oth- 


ers for multi-phase injection or to be used as the single-phase flooding liquid. 
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Downstream of the liquid selection valves there is a four-way valve that allows 
a Circulation loop through the core and another simultaneous loop back to a pump 
storage flask. The fractional flow compositions of the two loops are determined by 
the liquid selection valves and the flow rates specified by the pumps. Turning the 
four-way valve diverts the fluid from one loop to the other and vice versa. This 
permits an uninterrupted interchange of flow from the multi-phase injection stream 


to the single-phase flood stream at the start of the flood portion of the run. 


A pressure relief valve (set to 20 psi) along with a pressure gauge (0-200psi), 
points (6) and (7) on Figure 4.8, is installed just before the core/bypass inlet valve. 
The pressure relief valve is a safety measure to protect the pressure transmitter from 
possible damage due to over pressurization. The pressure gauge allows the system 
pressure to be read during pressure testing. Both of these devices have valves 
before them in order to isolate them from the system while the evacuated core is 


being flooded. This eliminates air from entering the core during this critical flood. 


All tubing used to connect the valves is either 1/8-inch stainless steel or 1/8- 
inch teflon tubing. The valves used are Whitey® two, three, or four-way, 1/8-inch 
swagelock® connecting valves. Between each pump outlet and the outlet control 
valve is a filter assembly and a pressure surge suppressor. A Q.5-micron sintered 
stainless-steel filter element is used in each filter assembly. Each surge suppressor 
consists of a vertical, graduated glass bulb attached to the outlet stream by a tee and 
1/4-inch adapter. Each bulb has an approximate volume of 6 mL. These suppres- 
sors are partially filled with fluid and the fluid levels recorded before and after each 


operation in the run. 


25 





26 


4.4 Porous Medium and Core Holder 
The porous medium consists of a strongly water-wet glass bead pack enclosed 
in a stainless steel core holder. The screen analysis of the number 8 glass beads 


used in run 2, and the number 9 glass beads used in runs 3 through 8 are shown in 


Table 4.2. 


Table 4.2 


Screen Analysis of Glass Beads 


mesh size | weight % no. 8 beads | weight % no. 9 beads 


0.54% <0.02% 
85.5% 72.8% 
13.81% 27.2% 
0.13% <0.02% 
0.02% 0.0% 





Two core holders were constructed of 1-inch O.D. stainless steel pipe, each 
36-inches long. One core is used in the apparatus at a time. The second core can 
be prepared for the next run during slow periods. Figure 4.9 shows the pressure tap 
and tubing configuration of the core. Holes for 1/8-inch stainless steel pressure taps 
were drilled symmetrically along the core. These pressure taps are silver soldered 
in place. 

The ends of each core are milled flat. One-inch swagelock® caps seal the 


ends of the core. Each cap had two holes drilled into it, one for the fluid stream 
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inlet or outlet and the other as a pressure port. The stainless steel tubing (1/8 inch) 
for the inlet or outlet duct and pressure port are silver soldered into these holes. A 
fine-mesh brass screen wrapped around a disc of 3/16-inch thick sintered brass mat- 
ting 1s placed in each core end cap between the end of the core and the core seat in 
the cap. This screen and disc keeps the glass beads in the core while allowing 


fluids to enter and leave the core with minimal pressure drop. 


Three-way valves located just before and just after the core allow the core to 
be bypassed. This allows the adjustment of fluid levels in the reservoir without 


having to pass these fluids through the core. 


Plexiglass® stands hold the bypass valves, core, and pressure tap valve mani- 


fold in order to give these portions of the system rigidity. 


4.5 Pressure Monitoring System 


Six symmetrically placed pressure taps along the core, along with unions and 
shut-off valves for each tap, make up the pressure tap manifold. Figure 4.9 shows 
the numbering of the pressure taps along the core. The pressure tap manifold is 
designed to allow the measurement of the pressure differential between an open 
pressure tap along the core and the end of the core. A Paroscientific Model 1020-D 
Digiquartz® Intelligent Transmitter is used to measure the pressure differential. 
The transmitter is connected to a microcomputer Бу a RS232 serial port card. The 
transmitter integrates the pressure readings over the sample time selected and 
corrects for any temperature changes during the sample period. The transmitter also 
has the capability to measure ambient temperature. Sampling speed (integration 
time) 1s set by adjusting the baud rate and pressure resolution. The baud rate for all 
experiments is set at 1200 and the pressure resolution is adjusted to change the 


sampling rate. For core linearity measurements and dynamic displacement 
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experiments the sample time is set at approximately one second. During steady- 


state runs, the sampling time is set at approximately thirty seconds. 


A BASIC program, (readpt.bas) contained in Appendix B, is used during 
experiments to instruct the transmitter to measure the initial ambient temperature, 
followed by continuous pressure readings. The program measures and records the 
elapsed time from when the program is initiated until the measurements are taken. 
When the program is given the stop command, the final temperature is measured. 
All pressure, temperature, and corresponding time measurements are displayed on 
the monitor and written to the microcomputer hard drive. A Paroscientific Model 
710 LCD visual display is used to display the pressure data generated by the Model 


1020-D transmitter. 


4.6 Collection System 


A three-way valve following the core/bypass outlet three-way valve, deter- 
mines if the fluid stream is directed to the fractional collector during dynamic dis- 


placement experiments or to the system reservoir during steady-state experiments. 


An ISCO Retriever II fraction collector is used to collect the fluids during 
dynamic displacement experiments. The Retriever II can collect fractions on the 
basis of time (0.1 to 999 minutes), counted drops (1 to 9990 drops) or other 
counted events (1-9990). For all runs, the time basis is used. The collection time 
is set to either 5 minutes or 3 minutes. Tubes used for collection in the Retriever II 
were 15 mL graduated centrifuge tubes. During a dynamic displacement when most 
of the core fluids have been displaced, and only small amounts of displaced fluids 
are being collected, the Retriever II collection method is replaced with hand 


changed 100 mL centrifuge collection tubes. 
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4.7 Additional Equipment 


In addition to the major equipment described so far, a few auxiliary items and 


pieces of equipment were required. 


In order to evacuate gases from the core prior to initial flooding, a vacuum line 
is connected by a tee between the core and the core/bypass outlet three-way valve 
(see figure 4.8). In series, between the core outlet tee and the vacuum pump, are a 
three-way valve, vacuum transducer, and 2000 mL vacuum flask. The three-way 
valve permits the vacuum pump to evacuate the core, isolate the core from the 
vacuum pump, or allow air into the evacuated flask. A Leybold Trivac "A" dual 
stage rotary vane pump, model D2A, is used as the vacuum pump. A Hastings 
vacuum gauge, model VT-6B, and vacuum tube, type DV-6R, is used to measure 


the vacuum drawn on the vacuum system. 


An IEC centrifuge, model HN-S II, is used to facilitate quicker separation of 


phases collected during the dynamic displacement experiments. 
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Chapter 5 


EXPERIMENTAL PROCEDURE 


5.1 System Core and Tubing Volume Determination 


Prior to any experimental runs, the volumes of tubings and cores are deter- 
mined. The core volumes are determined through linear measurements. The tubing 
volumes and core end cap volumes are determined by injecting distilled water into 
the tubing with a 3—mL syringe and reading the injected volume needed to fill the 


tubing. 


5.2 Pump Calibration 


Actual pump flow rates for each pump are determined by measuring the 
elapsed time, At, for a volume of fluid, V, to be displaced from the pump for a 
variety of specified flow rates. The actual flow rate for a particular specified flow 


rate 1S: 


V 
Чап = At (5.1) 


The specified flow rate versus actual flow rates for each pump are shown in Chapter 


6. 


5.3 Core Preparation 


The #8 glass beads for run 2 were sieved, and the 70-100 mesh beads used for 
the packing. All other runs used #9 glass beads that are washed in water, rinsed in 


distilled water and dried in a core drying oven for 48-72 hours at 50 ? c. 


The procedure used to pack the core is as follows: 
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1. The weights of the empty core holder and flask containing the glass beads are 


taken and recorded. 


2. The outlet end of the core 1s capped, and the glass beads poured into the core 
in 80 mL volume lifts. After each lift, a 3/4-in x 36-in aluminum rod is 
inserted into the core to hold the beads in place. The core is then tapped 
lightly up and down its length with a rubber headed hammer for several 


minutes in order to settle and compact the beads uniformly. 


3. The core is filled and tapped as described above until the beads are flush with 
the top of the core. An additional mound of beads approximately 1/8-inch in 
height is added to account for any additional settling. The inlet cap is then 


screwed into place and tightened. 


4. The weights of the full core holder and bead flask are again taken and 
recorded. After final tightening of the end caps, the threaded joint on each end 


cap, as well as each core/cap joint, is sealed with a fillet of vacuum grease. 


5.4 Determination of Porosity 


After each core is packed, and prior to conducting any experimental runs, the 


core's porosity is determined. By definition 


V 
pore 
= — 52 
ф vn. (5.2) 


The pore volume of each core glass bead pack is determined as follows: 


1. The core is first evacuated. This is accomplished by the following sequence of 
events. The core is isolated from the rest of the apparatus by closing all of the 
valves leading to the core. The valve leading from the core outlet to the 


vacuum pump is opened, and the core is then evacuated to a vacuum of not 
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greater than 100 microns of mercury. The vacuum pump normally runs over- 


night to achieve this level of vacuum. 


The apparatus is then set to flood the core with water while measuring the 
volume of water needed to accomplish the flood. This is accomplished by the 
following sequence of events. The core inlet and outlet valves are set to 
bypass. The liquid selection valves are set for a water flood. The water surge 
suppressor, pressure relief valve, and pressure gauge are isolated from the sys- 
tem by closing the appropriate valves. All other valves are then set to a no- 
flow position in order to minimize possible vacuum loss and/or fluid contami- 
nation of the core. The reservoir is filled from the water flask via the bypass 
with a minimum of one core bulk volume of water. This ensures that the injec- 
tion tubing network contains only water and that there is enough water in the 
reservoir to fill the core. The water pump feed valve is then set to draw water 
from the reservoir so that a loop is established from the reservoir to the pump 
through the bypass and back to the reservoir. The initial reservoir fluid levels 


and tubing heights, and surge suppressor fluid volumes are recorded. 


The core is then flooded. This is accomplished as follows: The water pump 
flow rate is set to 5.0 mL/min. The vacuum valve is closed and the core inlet 
valve is turned from bypass to the core, thereby pumping water from the reser- 
voir into the evacuated core. The pump is stopped at a volume several mL 
from the predicted volume needed to fill the core, and water is allowed to be 
pulled through the pump by the remaining vacuum. Once the water level has 
stabilized in the reservoir, the water surge valve is slowly and carefully opened 
to see if the reservoir water level drops indicating that a vacuum still remains. 


If the level drops, water is pumped one mL at a time into the core and the 
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5.5 


surge test repeated until there is no further drop in the surge fluid level. The 
core outlet valve is then opened and the system pressures and fluid levels 


allowed to equilibrate. 


The final reservoir and surge fluid levels, and tubing heights are recorded. 
Tubing displacement volume is determined by measuring the linear difference 
between the initial and final fluid heights multiplied by the number of tubes 
within the displaced zone and then multiplied by the cross-sectional area of the 
tubing. The pore volume is equal to the volume of the water that was 
pumped, minus the dead volume of the core under vacuum, minus the tubing 
displacement volume in the AH,O displaced, minus the change in the surge 


volume. 


AV (5.3) 


tubing — surge 


V pore = АУ во Б V dead N 


The porosity is then determined by the definition given in Equation (5.2). 


After the final fluid readings are taken, approximately 40 mL of water from the 
water flask is pumped through the core and into the fraction collector outlet to 
flush any fines from of the core. Unflushed fines could cause premature clog- 
ging of the system’s filters. The pressure tap manifold is then flushed and filled 


with water and the pressure transmitter connected. 


Determination of Core Linearity 


Horizontal steady-state linear flow through homogeneous porous media must 


have a linear pressure drop along the length of the porous media. Core pressure 


drop linearity is determined as follows: 


І 


Pressure transmitter sample rate is set to one sample per second. Water pump 


flow rate is set to less than the flow rate required to inject one pore volume 
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5.6 


per hour. The usual specified flow rate for linearity and absolute permeability 


determinations 1s 2.0 mL/min. 


The BASIC data recording program (readpt.bas) shown in Appendix B, is 
started and pressure drops are recorded across the core at each tap until the 
pressure stabilizes and a constant reading is observed. A normal recording 
time is 3-5 minutes per tap. The glass bead pack is used for experiments if 
the pressure drop across the core’s pressure taps are within 10% of each other. 
The pressure drops used for the test are across taps 2 and 3, AP» ,, taps 3 and 
4, АРа 4, and taps 4 and 5, AP, 5 (see figure 4.9). As only pressures across a 


specified tap and tap 6 can be measured, the equation 


AP AD Ape (5.4) 


a-b -~ 


is used to determine the pressure drop across the taps of interest. 


Determination of Absolute Permeability 


For horizontal steady-state linear flow it can be shown from Darcy’s Law that 


absolute permeability k,,, 1s 


] 
Kaps = Е (5.5) 


where; 


q: water flow rate in mL/sec 

p: water viscosity in ср at the final run temperature 
1; length of core in cm 

AP: final pressure drop across the entire core in atm 


A: cross-sectional area of core in cm? 
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5.7 Steady-State Experiments 


Steady-state experiments are run to achieve specific core saturations and deter- 
mine the effective and relative permeabilities at those saturations. A specific frac- 
tional fluid flow is injected into the core until the same fractional flow is produced 
from the core. Steady-state is assumed to be achieved at this point. Darcy’s law is 
used to determine the effective permeabilities and material balances are used to 


determine the core saturations. 


After the core characteristics are found, the core is brought to residual water 
saturation through a steady-state oil flood. From residual water saturation, the core 
is taken to a desired two-phase or three-phase saturation through a multi-phase 
steady-state flood. All steady-state experiments follow essentially the same pro- 
cedure, only the fractional flow of the injected fluids change from experiment to 


experiment. The following procedure is used for steady-state floods. 


1. The reservoir is charged, through the bypass, with enough of each fluid to 
reach the desired core saturation without expending all of any of the reservoir 


fluids being flooded. 


2. The apparatus valves are adjusted so that the flooding fluids are looped from 
the reservoir, to the pumps, through the multi-phase injection portion of the 


fluid injection valve system, on through the bypass and back to the reservoir. 


3. The pumps of the flooding fluids are set to achieve the desired fractional flow 
rates. The pumps are run until the multiphase injection, bypass, and reservoir 
return tubings are filled with the new fraction of fluids. All fluid levels in the 
reservoir and surge tubes, along with the tubing heights are taken and 


recorded. These readings are called the initial readings. 
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The pressure tap manifold valves are set to measure the pressure drop across 
the entire core (taps 1 and 6 are the only open valves). The pressure 
transmitter sampling rate is set. The rate used for all steady-state runs is one 
sample each thirty seconds. The BASIC data recording program, (readpt.bas) 
contained in appendix B, is begun and both core inlet and outlet valves are 


simultaneously turned from bypass to core. 


A minimum of three pore volumes of fluid is allowed to flow through the core. 
The pressure and reservoir fluid levels are monitored. Once the pressure and 
fluid levels have stabilized, steady-state is assumed to have been reached. The 
pumps are shut off and all valves switched to the no-flow positions. All fluid 
levels and tubing heights are again taken and recorded. These are known as the 


final readings. 


Saturations of the liquids and their effective and relative permeabilities are 


then determined. The following quantities are defined and used to determine satura- 


tions and permeabilities. 


|, 


Fractional flow of fluid i, f;; where q, and q; are the actual individual fluid 
flow rates and q, is the total flow rates of the pumps. 
shi Ч; 


Ё = ш ----- (5.6) 
| dt qi + qj 


Tubing volume, Viypings Where H; and H) are the top and bottom liquid inter- 
face heights in the reservoir measured in cm, Ay, is the cross-sectional area 
of the tubing, and Ny pes is the number of tubes, 


Viubing = (Hi — H3 ) Ambing Nubes (5.7) 


Initial and final fluid volumes; where Vyeading is the raw reservoir fluid column 


2" 





5.8 


volume reading in mL. 


MN Уор reading — Y bottom reading | V tubing (5.8) 


Change in surge volume; 


AV surge E V surge initial - V surge final (5:9) 


Change in fluid volume; 
A as па = У пыл = АУ, шее О Маса) (5.10) 


New core fluid volume, V; еу, соге; Where Vioidcore 18 the fluid volume in the 


core before the last flood, 


ее = М сос — AV; (5.11) 
Core saturations, S;; 
V. 
S. = 1lnew,core (5. 12) 
| V pore 


Effective permeability of each liquid at the final core saturation; 


qi Hy ! 
К = EET (5.13) 


Relative permeability of each liquid at the final core saturation; 


ki 
К; = k 





(5.14) 


abs 


Dynamic Displacement Experiments 


The dynamic displacement experiments are conducted to reduce the number of 


experiments needed to map the saturation region. In all dynamic displacement 


experiments a single phase is injected. Runs 2 and 3 end with a dynamic decane 


flood to verify the results with previous work. Runs 4 and 5 end with a dynamic 


water flood. In run 8, a dynamic decane flood is conducted prior to the dynamic 
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water flood in order to compare with the results of run 8’s earlier steady-state 


experiments. 
The following procedure is used for dynamic displacement experiments: 


1. Steady-state flow through the core is established to reach the desired initial 
core saturations. The injection into the core is through the four-way valve’s 
multiphase inlet and core injection outlet line. The valves for the single-phase 
injecting fluid are set to establish a loop from the flask, through the pump, into 
the four-way valve’s single-phase inlet and flask return outlet, and then back to 


the flask. 


2. A minimum of 5 pore volumes of flooding fluid is established in the flooding 
fluid flask. The specified flow rate of the flooding fluid pump is set to match 


the total specified flow rate of the steady-state flow rate. 


3. The pressure transmitter sampling rate is set to one sample per second. The 
BASIC data collection program (readpt.bas, found in Appendix B) is started 
and the fraction collector is set to the specified sample rate. In all runs except 
run 2, the fraction collector sample rate is set to 3 minutes. The 
reservoir/fraction collector selection valve is set to the fraction collector and at 
least 6 mL of the steady-state fluid fraction is pumped through the line to 
purge it. 

4. The dynamic flood is started by simultaneously turning the four-way valve to 
direct the single-phase flood to the core, advancing the fraction collector to the 
first dynamic displacement collection tube, and recording the elapsed time 


from the data collection program. 


5. The 15 mL collection tubes are used until the displaced core fluid volumes can 


no longer be read accurately. Normally this occurs when the oil volume is 





below 0.1 mL per tube. At this point the Retriever II fraction collector is shut 
off and 100 mL centrifuge tubes are placed by hand to collect the fluid. The 
empty and full weights of the 100 mL collection tubes and the elapsed collec- 
tion time of each tube are recorded. The filled fraction collection tubes are 
numbered and run through the centrifuge to separate the phases. The top fluid 
level reading for each fluid is recorded for each tube. The total fluid volume 
for the 100 mL tubes is determined by multiplying the weight difference of the 
full and empty tube by the density of the displacing fluid. 


Vi 100 mL ube = AWlibe Padiplacing fluid (5:15) 


The total fluid volume for the 100 mL tubes may also be determined by multi- 
plying the weight difference of the full and empty tube by the density of the 
displacing fluid if the flow rate during the collection is constant and equal to 


the calculated rate. 
V100 mL tube = Gacmal At (5.16) 


5.9 Fractional Recovery Analysis 


The raw fluid level data for each dynamic displacement experiment is used by 
a FORTRAN analysis program (datafrac.f), listed in Appendix B, to determine the 
actual fluid volumes in each tube, V;, the cumulative individual fluid collected 
volumes for each tube, V, the total cumulative fluid collected volume for each 
tube, V4, the cumulative pore volumes recovered for each fluid, R;, and the pore 


volumes injected , Q. Pore volumes recovered are defined as 


Уч 
R, = —— (5.17) 
V pore 


where V accounts for the fluid dead volumes in the core and tubing leading to the 
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fraction collector. 
n 
Va = DV; - fi Vaead (5.18) 
1 


Pore volumes injected are defined as 





Vu 
2) - (5.19) 
V pore 
where 
n 
Ma = УМ (5.20) 


The fluid recoveries are then plotted against the respective pore volume injected. 


The pressure data is converted by another FORTRAN program (pvst.f), con- 
tained in Appendix B, into the normalized pressure drop, ®, and plotted against the 


pore volumes injected, Q. The normalized pressure drop is defined as 


pues (5.21) 


Qcalc 





The program has the capability of averaging any specified number of the pressure 
data points in order to smooth the curve. Another FORTRAN program (variq.f), 
contained in appendix B, is a modified version of pvst.f that calculates the actual 
flow rates during the collection time of a collection vessel by dividing the actual 
collected volume by the time collected. This data is then used to determine the true 
normalized pressure drop and the actual pore volumes injected. This program is 


used if there is a varying flow rate during the dynamic flood experiment. 


4] 
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Chapter 6 


ERROR ANALYSIS 


In this section, the errors in calculating the following four quantities are dis- 
cussed: 


1. Porosity 
2. Saturation 
3. Absolute Permeability 


4. Relative Permeability 


The general approximation for the error is: 


_y |G] ay. 
AG = 2 m | AX; (6.1) 





where G = (хх) ХУ). 


6.1 Error Analysis for Porosity 


The equation for porosity calculation is: 


V 
DEC m (6.2) 
bulk 
where 
V pore = АУ но — Vaead 7 V mbing та AV surge (6.3) 


Vpore = VinH,O ^ VfinaH;O ^ Vdead 7 Vtubing ^ Vinsurge * V final,surge (6.4) 
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6 
= 2 
У рот = У НО Т V final, H,O = AV dead = Hi, E Нет Г, N= NS t ҰЛЫС тес 
i=] 


(6.5) 


where r, is the tubing radius and N is the number of tubes in the domain that is 


changing. 


Y bulk = Л po Есе (6.6) 
Пепсе: 


6 
2 
У; н.о — Vfna,H,O - 2,АУаеза — [Hana = Hin г М - Visurge + Vfinal,surge 
i=] 


К a, E 

(6.7) 
The summation term for the dead volume has six different values. Hence, the sum- 
mation is from 1 to 6. We will take the worst case of N = 4. We will compute 


the individual errors, EV апа ЕУ ло then we will compute the error in poros- 


pore? 
ity. 


6 
EV pore = EVinH,0 + EV final, H30 = 2 ЕД У dead 
i=] 


+ 47 E Ен, + EH na + Наш БЕ га 21, ` Er, + БУ шс + Еа ате 


(6.8) 
EVE = R 2 Icore (5544 Еее + Го El. (6.9) 
EV V EV EV 
Eo = E BV. = 6 I + = (6.10) 
Уык Мык У роге V bulk 
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The individual E values are controlled by the measurement techniques, and given 
below: 
Ето x0 5 cc 


ЕУ пано = 0.5 cc 
EAV gead = 0.02 cc 
ЕН, = 0.05 cm 
ЕНаа = 0.05 ст 
Ег, = 0.002 ст 
EVinsurge = 9.01 cc 
EV finalsurge = 0.01 cc 
Ercore = 0.002 cm 
Elore = 0.002 cm 


Typical measured values are: 


На + 0.0 cm 
Hina = 15.0 cm 
г + 0.3175 cm 

Toore = 1.0895 cm 


core 


] = 91.46 cm 


core 
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ф = 0.38 


Vouk = 341.0 cc 


134.0 cc 


< 
|| 


Applying the errors and typical values to Equations 6.8 and 6.9 yields: 


ЕУ. = 0.5 + 0.5 + 6: 0.02 
+ 4 n| 1.0895? 0.05 + 0.05 n 150 = 00 2. 1.0895 . 0.002] 


+ 0.01+ 0.01 = 1.968 = 1.97 сс 


po 


ЕУ, = л B 1.0895 - 91.46 - 0.002 4- 1.08952 0.002] 


= 1.259 = 1.26 cc 
Substituting EV.,. and EV}, into Equation 6.10 yields typical porosity error: 
97 1.26 
ENS US ee cm 
Фф | 134.0 = 
= 0.00699 


In percentage units = 0.7% error. 


6.2 Error Analysis for Saturation 


There are two different procedures for saturation determination: the steady- 


state and the non-steady-state. These two procedures are discussed next. 
6.2.1 Error Analysis for Steady-State Saturations 
The basic equation for average core saturations is: 


V: 
ны (6.11) 
роге 





o 


where Vj core 1S the volume of phase j in the core. Recall that EV,,,. = 1.97 cc. 


Now: 
Vi,core ы У оја, соге ре AV соге (6.12) 


Let us assume that for the first time fluid j is introduced to the core 


Wi old,core - 0 + 0.0. 


Therefore: 


Уса = - AV; (6.13) 


j,core 


АУ) Ра Міл — ‘j,final ~ "'jtubing - У ¡sin surge Е V j,final surge » f; Vaead (6.14) 
This is very similar to the Equation 6.4 except for the f; term. Let us take f; at the 


upper limit of 1. Then EV; = EV = 1.97 at the upper limit of complete 


j,core pore 


phase displacement. Hence: 


EV ј,соге S EV pore 
У; У 


j,core pore 








Typical saturation values are 0.1-0.8 and typical Vi, values are 


134.0: 0.1 – 134.0 · 0.8 = 13.4 - 107.2. 


Непсе, Гог S; = 0.1 


1.97 1.97 
. = 2 саа + Oooo amooo = А 1 2 SE T: 
ES; 0.1 um um 0.016 6% 
Еог S; = 0.8 
1.97 1.97 | 
. = 0$ | + —— | = 0.0264 = 2.6% 
Ерата p = р 


In summary, Е$; = 1%-3%. Рог са$е$ Шаг У; оасое * 0, the errors do not 


simply add up as the same fluid is used over and over. For example, the first oil 
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flood is done in a closed loop and Voiloldcore = 0. Now, the second steady-state 
oil-water flood (as in runs 2 and 3) is done without adding fluids to the cycling 
reservoir. Since the interfaces in the graduated cylinder are read in the same 
fashion every time, the errors cannot accumulate for one-phase independently. If 
the error increases for one phase it decreases for the other, hence, if the reading 
technique is constant, the errors should stay constant at about 1-2%. If new condi- 
tions are established in the cycling reservoir, the errors for phases existing in the 


core should double. 
6.2.20 Error Analysis for Non-Steady-State Saturations 


The determination of the saturation error for this case is rather difficult. The 
dynamic displacement saturation that we calculate using Ше Buckley-Leverett 


theory is assigned to the output end of the core and given by: 


ОК, 
5; = Sin F К, -Q 20 (6.16) 
| | OR; 
where R; is the recovery in pore volumes of phase j, and 3o = Rj. Let us 


assume that ES;;, is of the same order of the error we had before = 1-2%. The 
error in R; depends on the number of tubes we read and the individual error in the 
tubes. If we can reduce EV; pe = 0.1 cc then for the 15" tube (at about 
1.2 — 1.5 РУП), EV; total producea = 1.5 сс. This error can be reduced by knowing 


the overall injected volume. 


Since 


Yy. 
Да j,total produced (SOUS 


pore 
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ЕУ. , EV or 
ER, = R; СЕ Ее. (6.18) 
V pore V pore 
If a typical value for R; = 0.3 then: 
L5 1297 
ER; 2 0.3| ———— + —— | = 00156 = 1.6% 
j 0.3 - 134.0 rd i 
If the recovery is smaller, then the errors are smaller too. 
The term Q (PVD is: 
У; ; У 
Q = inj, total _ prod,total (6.19) 
У роте V pore 
and 
EV. . EV 
ns —— ADE ag (6.20) 
inj,total pore 








The error in total fluid produced can be reduced to = 0.5cc, at about 
] PVI ( = 150 cc), by considering pump rates and accumulations. At 1 PVI we use 
about 10-15 tubes. If we use these tubes for a PVI calculation, then 


EV prodtotal = 1.0 - 1.5 cc. We will use EV yoga ~ 0.5 cc. Hence: 


ot = 0.0184 = 1.8% 


EQ = 1.0 
Q ES 134.0 


rec 


PVI 





Typical slopes are on the order of 0.1 





| The error in the slope is difficult 
to assess as we are using eye construction. Let us assume that ER; = 0.01. 
Hence: 


ES; = (1% – 3%) + 1.0 % + (0.1): 18% +1:1% = 3.18% - 5.18% 





This error varies based upon what stage of the flood it is computed for, and can be 
reduced by consistent fluid level measurements and overall material balance correc- 


tions. 


6.3 Error Analysis for Absolute Permeability 


The error analysis for absolute permeability was documented by Obut (1989) 


and is summarized here. The absolute permeability is given by: 


ан! 
CS = = 72 
The error is defined by: 
Е Е. В EA EAP 
ЕК с = Kabs я + т + a + 7 + AP (0.23) 


The typical individual experimental measurements and their error values are: 


q = 2.00 cc/min 


Eq = 0.01 cc/min 


u = 1.00 cp 
Eu = 0.01 cp 
| = 90.00 cm 
El = 0.02 cm 
А = 3.729 ст? 


ЕА = 0.02 cm? 
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50 
АР = 1.0 ря 
EAP = 0.001 psi 


Kaps = 20 darcy 


_ 200 | 0:01 , 001, 0.02 , 002 , 0.01 


1 
Ек - 
a 2.0 1.0 900 3.729 10 





0.612 darcy: or = 3% error 


6.4 Error Analysis for Relative Permeability 


In this project, steady-state and dynamic relative permeabilities are calculated. 


We consider these two cases next. 


6.4.1 Error Analysis for Steady-State Relative Permeability 
The equation used for steady-state relative permeability of phase i is: 


К. 


1 





ki = в (6.24) 
abs 
The error is defined by: 
Ek. Ek 
Ек = s cp - (6.25) 
1 abs 


We assume that Ek, = Ek,p, as the k; equation is identical in format to the k,,, 


equation. If kj; = 0.1, Ek; = 0.0375 darcy, and k,,, = 20.0 darcy, then: 


0.612 0.612 
EEUU | 100 ' oe 


0.009 or = 9% 





6.4.2 Error Analysis for Dynamic Relative Permeability 


The Welge equation for the relative permeability as modified by Johnson, Bossler, 


and Naumann, (JBN) is: 

















К. 1 f: l M; 1 f: l u; 
Ko, = = а Т : J J = Е J J (6.26) 
abs abs A ] oP Kabs АФ 
а ох 
1 OP : | 
where a a> Ф 15 Ше distance between the intercept of the tangent with the 





y-axis, aS shown in Figure 2.1. Then: 


Ек Ef, ЕІ ЕШ EA E 
+ — + — + — + — + 
ro Е 1 А 





Ек, = ki == (6.27) 








Let us assume that E : Ф = 0.02 ®. The error in the fractional flow is the error 
in the slope of the recovery curve. We assumed before that this error is on the 
order of 1-296. The error in the viscosity is about 1%. The error in the area is 


about 0.1%. Hence: 


ЕК, = Қ; |00187 + 0.01 + 0.002 + 0.01 + 0.006 + 0.0? | 


К, ; [0.068 


= 6% 


The basic way to reduce this error is to reduce the error in the slopes. 


6.5 Relative Error Analysis 


The relative error does not consider constant errors that are added to the abso- 


lute error. For example, the error in measuring k,,, does not have any effect on the 


5] 





individual dynamic relative permeability errors of Equation 6.27. Only the actual 


measurements take part in the relative error. For example: 


E® 


Ef, 
Е.К; = 6, | + = (6.28) 
T 1 ‚1 | f; Ф 





Непсе: 


Е, К; = К; 10.03|= 3% 


These errors should be expected, and are reasonable in magnitude considering the 


equipment used for this project. 


6.6 Pump Flow Rates 


The error for the pump flow rates is discussed in Section 6.3. Figures 6.1 through 
6.3 are the relationships between the specified pump flow rates and the actual pump 
flow rates. The fluid flow path used in the determination of these relations is the 
following: flow from the fluid flask through the individual pump, through the core- 
bypass and out the fractional flow collection outlet. The dots in these figures 


represent the actual data points taken. 
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water pump actual flow rate [cc/min] 


Figure 6.1: 


= 
сл 


= 


> 
(л 





0 0.5 5 2 2.5 
water pump specified flow rate [cc/min] 


Water Pump - Actual vs. Specified Flow Rates 
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BA actual pump rate [cc/min] 


Figure 6.2: 
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decane pump actual flow rate [cc/min] 


Figure 6.3: 
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Chapter 7 


RESULTS 


Nine experimental runs were conducted in this research. A run consisted of a 
series of experiments using the same packed glass bead core. Run 1 was an equip- 
ment familiarity run. In run 1, the apparatus contained only distilled water. Multi- 
ple phase experiments began in run 2, consequently there are no multiphase results 


from run 1. 


Runs 2 and 3 were conducted to verify the apparatus. This was accomplished 
by running dynamic decane flooding experiments under conditions similar to those 
reported by Grader and O'Meara (1988). Run 2 was primarily a multiphase fami- 


liarity run. Run 3 was a repeat of Grader and O'Meara's (1988) run 15. 


Runs 4 and 5 were the first dynamic water flood runs. They were designed to 
enter the ternary diagram in the interior region of the irreducible water iso- 
saturation line. The objective was to provide insight into the general trends of 


dynamic water floods and give a background for selecting subsequent runs. 


Runs 6 and 7 were abandoned after initial water flooding because of non-linear 


pressure drops across the glass bead pack. 


Run 8 was designed to obtain data along the irreducible water iso-saturation 
line. The initial steady-state experiments were conducted to travel up and down the 
irreducible water iso-saturation line with the following objectives: to obtain a frac- 
tional fluid flow ratio versus core saturation plot, to examine the effects of hys- 
teresis, to provide benzyl alcohol and decane relative permeabilities, and to provide 


comparison with data obtained through dynamic displacement experiments. Run 8 
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concluded with a dynamic water flood originating at irreducible water saturation and 
residual oil saturation in the presence of decane. This provided data from a water 


flood at the upper limit of decane saturation. 


Run 9 was designed to obtain data along the two-phase water/benzyl alcohol 
saturation line (0% decane saturation), and along the residual decane iso-saturation 
line. These experiments provided data on the effect of decane in the system and 


provided data from a water flood at the lower limit of decane saturation. 


7.1 Three-Phase Decane Flood Experiments 


This section discusses the results of experimental runs 2 and 3. These two 
runs were used to test and fine-tune the apparatus. They were also used to validate 
the apparatus by comparing the results with previous work. Tables À.1 through A.4 
in Appendix À summarize the fluid recovery data and relative permeability data of 


these runs. 


Figure 7.1 is the pressure profile of the core linearity check for run 2. The 
plateaus in the profile show the pressure drop between the labeled core taps. In this 
linearity test each tap along the core was opened for 100 seconds. There is an 
unexplained discrepancy in this profile as the pressure drop between taps 2 and 6 
cannot be lower than the the drop between 3 and 6. The core is not completely 
linear as there is not equal drop between the three equally spaced taps, 2-3, 3-4, and 
4-5. The run was continued despite this nonlinearity just to gain familiarity with the 
equipment and to use the data as a comparison. There was no attempt to use this 


run for verification saturation mapping. 


Figures 7.2 and 7.3 show the recovery curves for the benzyl alcohol and the 


water with the injected decane. The initials of the recovered fluids identify the 
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recovery curves. The dots on the figures are the actual data points representing the 
data from each collection vessel. The data for these figures is found in Table A.1. 
Figure 7.2 shows the entire recovery curves for all three phases, while Figure 7.3 
enlarges the recovery curves of the water and alcohol. Breakthrough of the decane 
occurred at approximately 0.06 pore volumes injected as shown by the departure of 
the decane recovery curve from the x-axis at this point. The smooth transitions of 
both the water and alcohol curves after breakthrough indicate there is no fluid bank 
developed for either of these fluids during the flood and there will be little to no 
saturation shock jump. Why this condition developed is unclear and normally is 
related to fluid bypassing the glass bead pack along the interface of the case and 
glass bead pack. 


The saturation history of this run is shown in Figure 7.4. The glass bead pack 
was flooded with water (1), then flooded to residual water saturation with benzyl 
alcohol (2), and brought to a core saturation (3) of 49% water and 51% alcohol by 
a steady-state 1:1 water-benzyl alcohol flood. The saturation trajectory of the 
dynamic decane flood is shown by points (3) to (5). The Buckley-Leverett shock is 
shown by the jump from point (3) to point (4). The continuous saturation trajectory 
after breakthrough is shown by the line from point (4) to point (5). The dots along 
the trajectory are the saturations calculated at various pore volumes injected. The 
data for these points are contained in Table A.2 in Appendix A. The saturation tra- 
jectory was determined using the dynamic displacement trajectory technique 
presented in Chapter 2. 

The flow-rate normalized pressure drop versus the pore-volumes injected plot 
is shown as Figure 7.5. The pressure was normalized to the single calculated actual 


flow set on the decane pump. The sample rate was set at one sample per second 
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Figure 7.4: Saturation History and Trajectories; Run 2. 
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Figure 7.5: Dynamic Flood Normalized Pressure Profile; Run 2. 
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and one hundred data points were averaged to make one point on the graph. 


Figure 7.6 is the pressure profile of the core linearity check for run 3. The 
plateaus of the profile show the pressure drop between the labeled core taps. In this 
linearity test each tap along the core was opened for 5 minutes (300 seconds). 
Sample time was one second. The core had a maximum pressure difference of 
8.2% between zones AP3_, and AP, s. Other differences were 6.790 between AP» 4 


and АР, 4, апа 0.9% between АР, 3 апа АР, <. 


Figures 7.7 and 7.8 show the recovery curves for the benzyl alcohol and the 
water with the injected decane. The initials of the recovered fluids identify the 
recovery curves. The dots on the figures are the actual data points representing the 
data from each collection vessel. This data is contained in Table A.3 of Appendix 
A. Figure 7.7 shows the entire recovery curves for all three phases, while Figure 
7.8 enlarges the recovery curves of the water and alcohol. Breakthrough of the 
decane occurred at approximately 0.5 pore volumes injected as shown by the depar- 
ture of the decane recovery curve from the x-axis at this point. The abrupt decline 
of the slope of both the water and alcohol curves at breakthrough indicate there is 
no fluid bank developed for either of these fluids during the flood, however, there 


will be a significant Buckley-Leverett saturation shock. 


The saturation history of this run is shown in Figure 7.9. The glass bead pack 
was flooded with water (1), then flooded to residual water saturation, S,,,, with ben- 
zyl alcohol (2), and brought to a core saturation (3) of 63.6% water and 36.4% 
alcohol by a steady-state 10:1 water-alcohol flood. The saturation trajectory of the 
dynamic decane flood is shown by points (3) to (5). The Buckley-Leverett shock is 
shown by the jump from point (3) to point (4). The continuous saturation trajectory 


after breakthrough is shown by the line from point (4) to point (5). The dots along 
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Figure 7.9: Saturation History and Trajectories; Run 3. 





the trajectory are the saturations calculated at various pore volumes injected. This 


data is contained in Table A.4 of Appendix A. 


The flow-rate normalized pressure drop versus the pore-volumes injected plot 
is shown as Figure 7.10. The pressure was normalized to the single calculated 
actual flow set on the decane pump. The pressure sample rate was one sample per 
second. Every one-hundred pressure points were averaged and plotted to make the 


pressure profile. 


The water, oil and decane relative permeabilities for runs 2 and 3 are shown in 
semi-logarithmic plots in Figures 7.11 through 7.13. The dots on the curves are 
the actual data points calculated at various pore volumes injected. These data are 
contained in Tables A.2 and A.4 of Appendix A. Different sizes of glass beads 
were used between runs 2 and 3. Larger, size 8, glass beads were used for run 2. 
The resulting larger pore spaces and larger absolute permeability may account for 
the increase in oil and decane relative permeabilities and the resulting decrease in 


water relative permeability. 


The saturation trajectory and relative permeability trends of run 3 are con- 


sistent with the work of Grader and O’Meara (1988). 


7.2 Two-Phase Water-Benzyl Alcohol Relative Permeability Experiments 


The first portion of run 9 was dedicated to determining the two-phase relative 
permeabilities of a water-benzyl alcohol system. The run then goes on to determine 
the water-benzyl alcohol relative permeabilities at residual decane saturation. 
Tables A.5 through A.10 in Appendix A contain summaries of the data from this 


Tun. 


Figure 7.14 is the pressure profile of the core linearity check for run 9. The 
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Figure 7.10: Dynamic Flood Normalized Pressure Profile; Run 3. 
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Figure 7.11: BA Relative Permeabilities vs. BA Saturation ; Run 3. 
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Figure 7.12: Decane Relative Permeabilities vs. Decane Saturation; Run 3. 
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Figure 7.13: Water Relative Permeabilities vs. Water Saturation; Run 3. 
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plateaus of the profile show the pressure drop between the labeled core taps. In this 
linearity test each tap along the core was opened for five minutes (300 seconds). 
Sample time was one second. The core had a maximum pressure difference of 
4.6% between zones AP? 4 and AP, 5. Other differences were 0.95% between 
AP._3 and AP3_4, and 3.7% between AP4. 4 and АР, 5. 


The dynamic water displacement experiment to determine the water-benzyl 
alcohol relative permeabilities was run twice. The saturation history of the first 
experiment is shown in Figure 7.15. The glass bead pack was flooded with water, 
point (1), then flooded to residual water saturation, S,,, = 8.9 %, with the oil phase 
(benzyl alcohol) point (2). The oil saturated core was then dynamically water 
flooded to residual oil saturation (4), S,, = 11.6 % . Core-end oil saturation at 
breakthrough of the Buckley-Leverett shock was S, = 44.7 %, point (3). The dots 
along the trajectory are the saturations calculated at various pore volumes injected. 


These data are contained in Table A.5 of Appendix A. 


Figure 7.16 shows the recovery curve for the benzyl alcohol and the injected 
water. The initials of the recovered fluids identify the recovery curves. The dots 
on the figures are the actual data points representing the data from each collection 
vessel. These data are contained in Table A.6 of Appendix A. Breakthrough of the 
water occurred at approximately 0.625 pore volumes injected as shown by the 
departure of the water recovery curve from the x-axis at this point. The abrupt 
decline of the slope of the alcohol recovery curve at breakthrough indicates that no 
oil bank developed during the flood and accounts for the Buckley-Leverett satura- 


tion shock jump seen on the ternary saturation trajectory diagram, Figure 7.15. 


The pressure drop normalized to a fixed flow rate versus the pore volumes 


injected plot is shown as Figure 7.17. The pressure was normalized to the single 
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Figure 7.15: Saturation History and Trajectories; Run 9 - First Dynamic Water 
Flood 
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Figure 7.17: Normalized Pressure Drop (Fixed q) vs. PVI; Run 9 - First Dynamic Water Flood. 
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calculated actual flow set on the water pump. The pressure sample rate was one 
sample per second. Every one-hundred pressure points were averaged and plotted 
to make the pressure profile. A drop in the injection flow rate caused the sudden 
drop in pressure at approximately 1.2 pore volumes. This is believed to be due to a 
sudden sealing of the water supply flask stopper. This in turn hindered makeup air 
from entering the flask creating a vacuum and slowing the flow rate of the pump. 
The fixed flow rate used to calculate the normalized pressure drop could not com- 
pensate for this and therefore the dip in curve resulted. The FORTRAN program, 
pvst.f, used to calculate the ® versus Q data, was upgraded to correct this problem. 
This updated FORTRAN program (variq.f) contained in appendix B, calculates the 
actual flow rates for each collection period by dividing the the actual collected 
volume by the time collected. This data was then used to determine the true nor- 
malized pressure drop and the actual pore volumes injected. The results of this 
program's analysis is shown in Figure 7.18 and then plotted with the fixed-flow rate 
normalized pressure drop plot in Figure 7.19. As seen in Figure 7.19 the corrected 
plot is slightly above the fixed flow-rate plot, also, the corrected normalized pres- 


sure drop curve is much smoother than the uncorrected curve. 


The two-phase water-benzyl alcohol relative permeabilities of this first 
dynamic water flood are shown in Figure 7.20. The relative permeabilities were 
calculated using the variable flow-rate normalized pressure data. The dots on the 
curves are the actual data points calculated at various pore volumes injected. These 
data are contained in Table A.5 of Appendix A. 

A second dynamic water flood experiment of the oil saturated core at irreduci- 
ble water saturation was conducted to determine the effect of the flow-rate drop 


seen in the first dynamic water flood experiment. Figure 7.21 shows the saturation 
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Figure 7.19: Fixed vs. Variable Normalized Pressure Drop Comparison; Run 9 - First Dynamic Water Flood. 
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Figure 7.20: Oil and Water Relative Permeabilities vs. Water Saturation; 
Run 9 - First Dynamic Water Flood 
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Figure 7.21: Saturation History and Trajectories; Run 9 - Second Dynamic Water 
Flood 





history leading to the dynamic water flood and the saturation trajectory during the 
flood. The core was taken from residual oil saturation, point (1) Sj, = 8.45 %, at 
the end of the first water flood, to irreducible water saturation, point (2) 
S,, = 9.44 Фо, by an oil flood. The core was then dynamically water flooded with 
breakthrough, point (3), occurring at S, = 42.9 %. The flood ended at residual oil 
saturation, point (4) S,, = 9.16%. These data are contained in Table A.7 of 


Appendix A. 


Figure 7.22 shows the recovery curve for the benzyl alcohol and the injected 
water. The initials of the recovered fluids identify the recovery curves. The dots 
on the figures are the actual data points representing the data from each collection 
vessel. These data are contained in Table A.8 of Appendix A. Breakthrough of the 
water occurred at approximately 0.625 pore volumes injected as shown by the 
departure of the water recovery curve from the x-axis at this point. The abrupt 
decline of the slope of the alcohol recovery curve at breakthrough indicates that no 
oil bank developed during the flood and accounts for the Buckley-Leverett satura- 


tion shock jump seen on the ternary saturation trajectory diagram (Figure 7.21). 


The pressure drop normalized to a fixed flow rate versus the pore volumes 
injected plot is shown as Figure 7.23. The pressure was normalized to the single 
calculated actual flow rate set on the water pump. The pressure sample rate was 
one sample per second. Every one-hundred pressure points were averaged and plot- 
ted to make the pressure profile. The normalized pressure was recalculated using 
the FORTRAN program variq.f (contained in Appendix B), and the data plotted in 
Figure 7.24. 


The two-phase water-benzyl alcohol relative permeabilities of this second 


dynamic water flood are shown in Figure 7.25. The relative permeabilities were 
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Figure 7.23: Normalized Pressure Drop (Fixed q) vs. PVI; Run 9 - Second Dynamic Water Flood. 
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Figure 7.24: Normalized Pressure Drop (Variable q) vs. PVI; Run 9 - Second Dynamic Water Flood. 
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Figure 7.25: Oil and Water Relative Permeabilities vs. Water Saturation 
Run 9 - Second Dynamic Water Flood. 





calculated using the variable flow rate normalized pressure data shown in Figure 
7.24. The dots on the curves are the actual data points calculated at various pore 
volumes injected. These data are contained in Table A.7 of Appendix A. A com- 
parison of run 9’s two-phase water-oil relative permeability data from the first and 
second dynamic water flood is shown in 7.26. Curves A and D are the second 
water flood’s relative oil permeabilities and relative water permeabilities respec- 
tively, while curves B and C are the first water flood’s relative oil permeabilities 


` and relative water permeabilities respectively. 


The last phase of run 9 consisted of a third dynamic water flood but at residual 
decane saturation. Figure 7.27 shows the saturation history leading to this third 
dynamic water flood and the saturation trajectory during the flood. The core was 
taken from 100% water saturation, point (1) prior to flood 1, to irreducible water 
saturation, point (2) St, = 9.44 %, through the series of water and oil floods out- 
lined in the description of dynamic water floods 1 and 2. The glass bead pack was 
then brought to point (3) by a steady-state decane flood, Sy = 74.8 %, and then 
back to point (4) by a steady-state oil flood. Point (4) represents residual water, 
Sw = 9.5%, and residual decane, Ха + 8.1%, in the presence of oil, 


S, = 82.4 %. 


The third dynamic water flood was conducted with breakthrough (5) occurring 
at S, = 33.8 %, and a continuous saturation trajectory at residual decane satura- 
tion from point (5) to point (6), 5, - 83.6 % and Sœ = 8.29 %. Decane was not 
produced during this flood. These data are contained in Table A.9 of Appendix A. 

Figure 7.28 shows the recovery curve for the benzyl alcohol and the injected 


water. The initials of the recovered fluids identify the recovery curves. The dots 


on the figures are the actual data points representing the data from each collection 
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Figure 7.26: Relative Permeability Comparison; 
Run 9 - First and Second Dynamic Water Floods. 
A = Oil Relative Permeability, Second Dynamic Flood 
B = Oil Relative Permeability, First Dynamic Flood 
C = Water Relative Permeability, First Dynamic Flood 
D = Water Relative Permeability, Second Dynamic Flood 
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Figure 7.27: Saturation History and Trajectories; Run 9 - Third Dynamic Water 
Flood 
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vessel. These data are contained in Table A.10 of Appendix A. Breakthrough of 
the water occurred at approximately 0.54 pore volumes injected as shown by the 
departure of the water recovery curve from the x-axis at this point. The abrupt and 
continuous decline of the slope of the alcohol recovery curve at breakthrough indi- 
cates that no oil bank developed during the flood and accounts for the Buckley- 
Leverett saturation shock jump seen on the ternary saturation trajectory diagram, 


Figure 7.27. 


The plot of pressure drop normalized to a fixed flow rate versus the pore 
volumes injected is shown as Figure 7.29. The pressure was normalized to the sin- 
gle calculated actual flow set on the water pump. The pressure sample rate was one 
sample per second. Every one-hundred pressure points were averaged and plotted 
to make the pressure profile. The normalized pressure was recalculated with the 
actual varying injected flow rate using the FORTRAN program (variq.f) contained 
in Appendix B. This data is shown in Figure 7.30. A comparison of the two cal- 
culated normalized flow rates is shown in Figure 7.31. The agreement of the pres- 
sure profiles is much better than the comparison of pressure profiles in the first 
water flood shown in Figure 7.19. This is because the flow rates in the third water 
flood are more consistent and closer to the calculated fixed flow rate of the water 
pump. This shows that the fixed flow rate analysis used in subsequent data presen- 
tations is valid as long as there is no significant deviation from the predicted pump 


flow rate. 


The two-phase water-benzyl alcohol relative permeabilities measured at resi- 
dual decane saturation during the third dynamic water flood are shown in Figure 
7.32. The relative permeabilities were calculated using the variable flow rate nor- 


malized pressure data shown in Figure 7.30. The dots on the curves are the actual 
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Figure 7.31: Fixed vs. Variable Normalized Pressure Drop Comparison; Run 9 - Third Dynamic Water Flood. 
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Figure 7.32: Oil and Water Relative Permeabilities vs. Water Saturation; 
Run 9 - Third Dynamic Water Flood. 





data points calculated at various pore volumes injected. These data are contained in 
Table A.9 of Appendix A. A comparison of run 9’s two-phase water-oil relative 
permeability data from the second dynamic water flood and from the third water 
flood at residual decane saturation is shown in Figure 7.33. Curves A and C are the 
second water flood’s oil relative permeabilities and water relative permeabilities, 
respectively, while curves B and D are the third water flood’s oil relative per- 
meabilities and water relative permeabilities, respectively. As expected, relative 
permeabilities of the oil and water are lower in the third water flood because of the 
presence of residual decane. Note also, the oil relative permeability is reduced 
more significantly then the water relative permeability by the presence of residual 


decane Saturation. 


7.3. Two-Phase Benzyl Alcohol-Decane Relative Permeability Experiments 


The first portion of run 8 was designed to determine the two-phase relative 
permeabilities of a benzyl alcohol-decane system at residual water saturation. This 
was accomplished through a series of steady-state experiments and a dynamic 
decane flood. The run concluded with a dynamic water flood of the glass bead 
pack at residual water saturation and residual oil saturation in the presence of 


decane. 


Figure 7.34 is the pressure profile of the core linearity check for run 8. The 
plateaus of the profile show the pressure drop between the labeled core taps. In this 
linearity test each tap along the core was opened for five minutes (300 seconds). 
Sample time was one second. The core had a maximum pressure difference of 
9.406 between zones АР; ҙ апа АР, 5. Other differences were 3.096 between AP» 4 


and АРа 4, and 6.6% between AP3_4 and AP, >. 
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Figure 7.33: Relative Permeability Comparison; 
Run 9 - Second and Third Dynamic Water Floods. 
A = Oil Relative Permeability, Second Dynamic Flood 
B = Oil Relative Permeability, Third Dynamic Flood 
C = Water Relative Permeability, Second Dynamic Flood 
D = Water Relative Permeability, Third Dynamic Flood 
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The overall saturation path of run 8 is shown in Figure 7.35. The glass bead 
pack was flooded with water, point (1), then flooded to residual water saturation 
with the oil phase (benzyl alcohol), point (2). The oil saturated core was then 
brought to point (3) by a series of steady-state decane/oil floods, called the drainage 
series. Point (3) represents residual water saturation and residual oil saturation in 
the presence of decane. The glass bead pack was then brought down to residual 
water saturation and residual decane saturation in the presence of oil, point (4), by a 
series of decane/oil floods in the reverse sequence of fractional flow rates that made 
up the drainage series of experiments. The series of experiments from point (3) to 
point (4) is called the imbibition series. To complete the run the glass bead pack 
was brought back to residual water saturation and residual oil saturation in the pres- 


ence of decane, point (5), and then dynamically water flooded to point (6). 


Figure 7.36 shows the saturation history of the steady-state drainage series of 
experiments from run 8. Table 7.1 lists the saturations of each point, the fractional 
flow ratios and actual specified flows that were used to reach the point, and the 
relative permeabilities at the point. Figure 7.37 plots the fractional flow of oil and 
the resulting glass bead pack saturations obtained from the results of these steady- 
state drainage series experiments at residual water saturation. The solid dots are the 


actual calculated data points. These data are contained in Table 7.1 . 


The relative permeability data for oil and decane versus oil saturation from 
these drainage series of experiments is plotted in Figure 7.38. The solid dots are 
the actual calculated data points from each experiment and are contained in Table 


Tels 


Figure 7.39 shows the saturation history of the steady-state imbibition series of 


steady-state experiments from run 8. Table 7.2 lists the saturations of each point, 
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Figure 7.35: Overall Saturation History; Run 8. 
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Figure 7.36: Steady-State Drainage Saturation History; Run 8. 
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Table 7.1 


Summary of Run 8 Drainage Experiments 


dec · oil 
[mL/min] 
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Figure 7.39: Steady-State Imbibition Saturation History; Run 8. 
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Table 7.2 


Summary of Run 8 Imbibition Experiments 
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the fractional flow ratios and specified flow rates that were used to reach the point, 
and the relative permeabilities at the point. The relative permeability data for oil 
and decane versus oil saturation from these imbibition series of experiments is plot- 
ted in Figure 7.40. The solid dots are the actual calculated data points from each 
experiment. The comparison of the drainage (Small dots) and imbibition (large 
dots) relative permeabilities is shown in Figure 7.41. Figure 7.41 shows relatively 
unchanged values for benzyl alcohol relative permeabilities, while the decane curves 


show some hysteresis effects by shifting their slopes. 


The run 8 two-phase decane-benzyl alcohol experiments at residual water 
Saturation concluded with a decane flood of the glass bead pack at residual water 
saturation, residual decane saturation in the presence of oil. The saturation trajec- 
tory of this experiment is projected in Figure 7.42. The starting point of the 
dynamic decane flood is point (1) representing residual water saturation, 
S;w = 3.1 %, and residual decane saturation, S,; = 18.1 %, in the presence of oil. 
The Buckley-Leverett shock breakthrough occurred at point (2), Sy = 22.1 %, the 
saturation trajectory continues until point (3) representing residual water saturation, 
Sw = 3.1 %, and residual oil saturation, S, = 14.7 %, in the presence of decane. 
The dots are the actual calculated data points from the decane flood experiment and 


are contained in Table A.11 contained in Appendix A . 


Figures 7.43 and 7.44 show the recovery curves for the benzyl alcohol with 
the injected decane. The initials of the recovered fluids identify the recovery 
curves. The dots on the figures are the actual data points representing the data from 
each collection vessel. These data are contained in Table A.12 of Appendix A. 
Figure 7.43 shows the entire recovery curves for both phases, while Figure 7.44 


enlarges the recovery curve of the benzyl alcohol. Breakthrough of the decane 
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Figure 7.40:  Decane and Oil Relative Permeabilities vs. Oil Saturation; 
Run 8 - Steady-State Imbibition Series. 
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Run 8 - Drainage and Imbibition Series. 
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Figure 7.42: Saturation History and Trajectories; 
Run 8 - Dynamic Decane Flood. 
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occurred at approximately 0.02 pore volumes injected as shown by the departure of 
the decane recovery curve from the x-axis at this point. The smooth transition of 
the benzyl alcohol curve after breakthrough indicates that no oil bank developed 
during the flood and accounts for the small saturation shock jump seen in moving 
from point (1) to point (2) on Figure 7.42. 

The plot of pressure drop normalized to a fixed flow rate versus the pore 
volumes injected is shown as Figure 7.45. The pressure was normalized to the sin- 
gle calculated actual flow set on the decane pump. The pressure sample rate was 
one sample per second. Every one-hundred pressure points were averaged and plot- 
ted to make the pressure profile. Figure 7.46 is an enlargement of the early portion 
of Figure 7.45. As seen in Figures 7.43 and 7.44, breakthrough occurred very 
early at 0.02 pore volumes injected. Figure 7.46 makes it easier to evaluate the 
slope of the normalized pressure drop curve when determining relative permeabili- 


ties. 


The two-phase decane-benzyl alcohol relative permeabilities measured at resi- 
dual water saturation obtained from the drainage (small dots) and imbibition (large 
dots) series of steady-state experiments and the dynamic decane flood are shown in 
Figures 7.47 and 7.48. The steady-state data points are the stand alone points 
while the relative permeability points calculated from the dynamic flood data are 
connected by a curve. The dynamic flood relative permeabilities were calculated 
using the fixed flow rate normalized pressure data shown in Figures 7.45 and 7.46. 
The dots on the curves are the actual data points calculated at various pore volumes 
injected. These data are contained in Table 7.1 and 7.2 in this chapter and A.11 of 
Appendix A. The oil relative permeabilities are approximately the same in both the 


steady-state and dynamic floods, where as the decane relative permeabilities are 
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Figure 7.47: Comparison of Dynamic and Steady-State Calculated Relative Decane 
Permeabilities; Run 8. 
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Figure 7.48: Comparison of Dynamic and Steady-State Calculated Relative Oil 
Permeabilities; Run 8. 
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consistently higher in the steady-state experiments. The reason for the disparity in 


the decane relative permeabilities is unclear and requires more research. 


7.4 Three-Phase Water Flood Experiments 


This section discusses the results of the three-phase dynamic water flood 
experiments. These are conducted in runs 4 and 5, and at the conclusion of run 8. 
This series of runs is designed to map the interior region of the ternary diagram. 
This is attempted by making a series of water floods originating from several points 
spaced along the residual water saturation line at varying decane and benzyl alcohol 


(oil) saturations. 


Figure 7.49 is the pressure profile of the core linearity check for run 4. The 
plateaus of the profile show the pressure drop between the labeled core taps. In this 
linearity test each tap along the core was opened for 5 minutes (300 seconds). 
Sample time was one second. The core had a maximum pressure difference of 
7.290 between zones AP» з and AP, 5. Other differences were 1.670 between AP» 4 


апа АР; 4, and 5.796 between AP4. 4 and AP, 5. 


Figures 7.50 and 7.51 are the run 4 recovery curves for benzyl alcohol and 
water with the injected decane. The initials of the fluids identify the recovery 
curves. The dots on the figures are the actual data points representing the data from 
each collection vessel. These data are contained in Table A.13 of Appendix A. 
Breakthrough of the water occurred at approximately 0.39 pore volumes injected as 
shown by the departure of the water recovery curve from the x-axis and the abrupt 
changes in the slopes of the decane and benzyl alcohol recovery curves. The almost 
vertical increase of the slope of the alcohol recovery curve at breakthrough indicates 


that an oil bank developed during the flood. The tangent to the recovery curve has 
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a negative y-intercept which indicates an increase in core-end oil saturation at 
breakthrough. The zero slope of the decane recovery curve after breakthrough indi- 
cates that the decane was removed from the core with piston like displacement 
down to a residual decane saturation prior to breakthrough and, therefore, no decane 


was produced after breakthrough. 


The saturation history of run 4 is shown in Figure 7.52. The glass bead pack 
was flooded with water, point (1), then flooded to residual water saturation, 
S,, 7 0.99 95, with the oil phase (benzyl alcohol) at point (2). The core at residual 
water saturation was flooded in a steady-state experiment using a specified frac- 
tional flow of 10:1 decane to oil to point (3). The saturations at point (3) are 
Sa = 495%, S, = 49.5 %, and S,, = 1.0 %. The glass bead pack was then 
dynamically water flooded. Core-end oil saturation at breakthrough of the Buckley- 
Leverett shock was S, = 85.0 %, point (4). The line from point (4) to point (5) 
represents the continuous saturation trajectory between breakthrough and the final 
flood conditions at residual oil and decane saturations. The points are the satura- 
tions calculated at various pore volumes injected. These data are contained in 
Table A.14 of Appendix A. The jump in oil saturation values at breakthrough as 
predicted by the recovery curve of Figure 7.51 is clearly obvious in the oil satura- 
tion jump from point (3) to point (4) in Figure 7.52. 

The pressure drop normalized to a fixed flow rate versus the pore volumes 
injected plot for run 4’s dynamic water flood is shown as Figure 7.53. The pressure 
was normalized to the single calculated actual flow set on the water pump. The 
pressure sample rate was one sample per second. Every one-hundred pressure 


points were averaged and plotted to make the pressure profile. 


Figure 7.54 is the pressure profile of the core linearity check for run 5. The 
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Figure 7.52: Saturation History and Trajectories; Run 4. 
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plateaus of the profile show the pressure drop between the labeled core taps. In this 
linearity test each tap along the core was opened for 5 minutes (300 seconds). 
Sample time was one second. The core had a maximum pressure difference of 
4.2% between zones AP» 4 and AP, s. Other differences were 3.996 between AP» 4 


апа АР; 4, and 0.4% between АР; 4 and AP, ;. 


Figures 7.55 and 7.56 are the run 5 recovery curves for benzyl alcohol and 
decane with the injected water. The initials of the fluids identify the recovery 
curves. The dots on the figures are the actual data points representing the data from 
each collection vessel. These data are contained in Table A.15 of Appendix A. 
Breakthrough of the water occurred at approximately 0.46 pore volumes injected as 
shown by the departure of the water recovery curve from the x-axis and the very 
abrupt changes in the slopes of the decane and benzyl alcohol recovery curves. The 
almost vertical increase of the slope of the alcohol recovery curve in Figure 7.56 at 
breakthrough indicates that an oil bank developed during the flood. The tangent to 
the recovery curve has a negative y-intercept which indicates an increase in core- 
end oil saturation at breakthrough. The zero slope of the decane recovery curve 
after breakthrough indicates that the decane was removed from the core with piston 
like displacement down to a residual decane saturation prior to breakthrough and 


hence, none was produced after breakthrough. 


The saturation history of run 5 is shown in Figure 7.57. The glass bead pack 
was flooded with water, point (1), then flooded to residual water saturation, 
Siw = 2.1 %, with the oil phase (benzyl alcohol) at point (2). The core at residual 
water saturation was steady-state flooded with a specified fractional flow of 20:1 
decane to oil to point (3). The saturations at point (3) are Sy = 58.5 %, 


5, = 39.4 %, and S,,, = 2.1 %. The glass bead pack was then dynamically water 





129 


`$ UNY залтмо) Алалод9ч рооа опишеш а :66'/ ат 4 


(разоэ[и! эшпол элоа) © 
ӯ $ © 






(ла) Алелозе: 


со 





130 


"с ий ‘SIANO ÁJIAOIDA POO|4 општа рориваха :9С'/ 211814 


(рәјоә[и эшптол элоа) © 
v € © 





ч 
o 


50 


(Ла) Аләлоэәл 





131 


ОЕСАМЕ 
100 


100 100 
МАТЕВ 50 2 OIL 


Figure 7.57: Saturation History and Trajectories; Run 5. 
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flooded. Core-end oil saturation at breakthrough of the Buckley-Leverett shock was 
So = 82.8 %, point (4). The line from point (4) to point (5) represents the con- 
tinuous saturation trajectory between breakthrough and the final flood conditions at 
residual oil and decane saturations. The dots are the saturations calculated at vari- 
ous pore volumes injected. These data are contained in Table A.16 of Appendix A. 
The oil bank saturation jump predicted by the recovery curve of Figure 7.56 is 
clearly obvious in the oil saturation jump from point (3) to point (4) in Figure 7.57. 
This trajectory, and that in Figure 7.52, differ in two ways from the water flood of 
the two-phase water-oil experiment at residual decane saturation in run 9, shown in 
Figure 7.27. The first difference is that there is a Buckley-Leverett jump along the 
constant decane Saturation line in Figure 7.27, and the Buckley-Leverett shock in 
Figures 7.52 and 7.57 occurs along the residual water saturation line. This means 
that there is a continuous saturation trajectory along the residual decane saturation 
line in Figures 7.52 and 7.57. Along this line there can be continuous relative per- 
meability measurements. In Figure 7.27 the Buckley-Leverett jump prohibits a con- 
tinuous measurement of relative permeabilities. The second difference is that the 
residual decane saturation line in Figure 7.27 is lower than the residual decane 


saturation lines in Figures 7.52 and 7.57. 


The pressure drop normalized to a fixed flow rate versus the pore volumes 
injected plot for run 5’s dynamic water flood is shown as Figure 7.58. The pressure 
was normalized to the single calculated actual flow set on the water pump. The 
pressure sample rate was one sample per second. Every one-hundred pressure 


points were averaged and plotted to make the pressure profile. 


Oil relative permeabilities calculated from the data of the dynamic water floods 


of runs 4 and 5 are plotted in Figure 7.59. Water relative permeabilities calculated 
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Figure 7.58: Normalized Pressure Drop vs. PVI; Run 5. 
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Figure 7.59: Oil Relative Permeabilities; Runs 4 and 5. 
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from the data of the dynamic water floods of runs 4 and 5 are plotted in Figure 
7.60. As can be seen from these two figures, there is a very good match of the 
permeability data from runs 4 and 5. Figure 7.61 shows the oil and water relative 
permeabilities from the dynamic water floods of runs 4 and 5 in the conventional 
Cartesian plot of relative permeabilities. The points in Figures 7.59 through 7.61 
are the actual relative permeability values obtained from the data at various pore 
volumes injected. The data for Figures 7.59 through 7.61 is listed in Table A.14 
and A.16 of Appendix A. 


As discussed in Section 7.3, the final portion of run 8 was a dynamic water 
flood of the glass bead pack from the initial conditions of residual water and oil 
Saturations in the presence of decane down to the final conditions of residual oil and 
decane in the presence of water. The initial and final points of this dynamic water 


flood where shown to be points (5) and (6) respectively on Figure 7.35. 


Figures 7.62 and 7.63 are the run 8 dynamic water flood recovery curves for 
benzyl alcohol and decane with the injected water. The initials of the fluids iden- 
tify the recovery curves. The dots on the figures are the actual data points 
representing the data from each collection vessel. These data are contained in 
Table A.17 of Appendix A. Breakthrough of the water occurred at approximately 
0.56 pore volumes injected as shown by the departure of the water recovery curve 
from the x-axis and the very abrupt change in the slope of the decane recovery 
curve. The small increase of the slope of the alcohol recovery curve at breakthrough 
indicates that a small oil bank developed during the flood. The tangent to the 
recovery curve has a negative y-intercept which indicates a slight increase in core- 
end oil saturation at breakthrough. As with the recovery curves of the decane in 


runs 4 and 5, the zero slope of the decane recovery curve after breakthrough indi- 
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Figure 7.60: Water Relative Permeabilities; Runs 4 and 5. 
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cates that the decane was removed from the core with piston like displacement 
down to a residual decane saturation prior to breakthrough, and therefore, no decane 
was produced after breakthrough. It should be noted that the oil bank formed in 


runs 4 and 5 is more significant than the bank of run 8. 


The saturation history of run 8’s dynamic water flood is shown in Figure 7.64. 
The glass bead pack was taken from the final conditions of residual oil and water in 
the presence of decane after the run 8 dynamic decane flood, point (1), to residual 
oil and decane saturation in the presence of water, point (3) through a dynamic 
water flood. Core-end oil saturation at breakthrough of the Buckley-Leverett shock, 
point (2), was S, = 20.7 %. The line from point (2) to point (3) represents the 
continuous saturation trajectory between breakthrough and the final flood conditions 
at residual oil and decane saturations. The small oil bank saturation increase 
predicted by the recovery curve of Figure 7.63 is seen in the oil saturation jump 
from point (1) to point (2). No decane was produced between points (2) and (3). 
The dots represent the saturations calculated at various pore volumes injected. 


These data are contained in Table A.18 of Appendix A. 


The pressure drop normalized to a fixed flow rate versus the pore volumes 
injected plot for run 8’s dynamic water flood is shown as Figure 7.65. The pressure 
was normalized to the single calculated actual flow set on the water pump. The 
pressure sample rate was one sample per second. Every one-hundred pressure 
points were averaged and plotted to make the pressure profile. The very small 
slope of the curve after breakthrough corresponds to the very little amount of oil 


produced after breakthrough. 
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Figure 7.64: Saturation History and Trajectories; Run 8 - Dynamic Water Flood. 
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7.5 Unacceptable Core Linearity Runs 


Runs 6 and 7 were not continued after initial core flooding because of core 
linearity problems. Figure 7.66 is the pressure profile of the core linearity check for 
run 6. The plateaus of the profile show the pressure drop between the labeled core 
taps. In this linearity test each tap along the core was opened for five minutes (300 
seconds). Sample time was one second. The core had a maximum pressure 
difference of over 40096 between zones AP, 4 and AP, 5. This exceeded the 
allowed 10% difference. 


Figure 7.67 is the pressure profile of the core linearity check for run 7. The 
plateaus of the profile show the pressure drop between the labeled core taps. In this 
linearity test each tap along the core was opened for five minutes (300 seconds). 
Sample time was one second. The core had a maximum pressure difference of over 


220% between zones AP._3 and AP,_5. This exceeded the allowed 10% difference. 


The permeability nonlinearity of these packs is probably caused by settling of 
the glass beads and fluid bypassing the glass bead pack along the pack/case inter- 
face. This explanation is further substantiated by the following two observations. 
While flowing water through the core, tapping the core lightly produced a reduction 
in the pressure differential read out, this could be caused by unsettling the glass 
bead pack. When the cores were opened to remove the glass beads and clean the 
core, there was a void approximately 1/4-inch deep at the upstream end of the core, 


this would have allowed the glass bead packs to expand and unsettle. 
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Chapter 8 


DISCUSSION 


8.1 Verification 


The initial portion of the research is concerned with decane flooding and the 
verification of previous work. This work is accomplished in runs 2 and 3. The 
uncertainty in the linearity of the core used in run 2 (see Figure 7.1), resulted in run 
3 being used to verify the apparatus. Run 3 reproduced results and trends of 
Grader and O’Meara (1988). The initial core saturation of run 3 is very close to 
run 15 of the work reported by Grader and O’Meara (1988). The saturation trajec- 
tory of run 3 (Figure 7.9) indicates a good match with the saturation trajectory of 
their (1988) run 15. This comparison of saturation trajectories is shown in Figure 
8.1. In Figure 8.1, the saturation trajectories of run 3 and Grader and O’Meara’s 


run 15 are represented by curves A and B respectively. 


8.2 Apparatus 


The apparatus is capable of the flexibility that was designed into it. Both 
steady-state and non-steady-state experiments are conducted under a variety of 
injected fluid scenarios. The apparatus showed the ability to reproduce results. An 
example of this is seen by the repeat of run 9's first dynamic water flood in the 
second run 9 water flood. The comparison of the relative permeabilities of these 
two runs (Figure 7.26), the comparison of the pressure drops normalized to the vari- 
able flow rate (Figures 7.18 and 7.24), and the comparison of the recovery curves 


(Figures 7.16 and 7.22) indicate close agreement. 
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Figure 8.1: | Comparison of Saturation Trajectories; 
Run 3 (marked as A) vs. Grader and O’Meara (1988) Run 15 
(marked as B). 
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Some limitations shown by the apparatus are seen in the inability to obtain a 
continuous three-phase recovery curve. This differs from the pressure recovery 
curve where a data point 15 obtainable every second with very good resolution. The 
pressure recovery curve is nearly continuous. A source of error and a bottleneck in 
the data analysis is the fact that the data analysis is not automated. Recovery curve 
and normalized pressure curve slopes are taken by hand. Future improvements to 
the apparatus should include the addition of these modifications to expedite the 


analysis and also reduce some of the error associated with hand analysis. 


Achieving uniform packing of the glass bead core is the most difficult step in 
the experimental procedure. The acceptable linear pressure drop of an individual 
core did not mean that the core is the same in terms of porosity and absolute per- 
meability as other cores. This уапапсе among the cores should account for most of 
the difference between similar terms calculated in different runs. It cannot account 
for the difference of similar terms calculated within the same run, 1.e. using the 
same core. Table A.19 of Appendix A lists some of the core characteristics of each 


run. 


8.3 Exploratory Findings 


The dynamic water floods of runs 4, 5, 8, and 9 are the new areas of investiga- 
tion conducted by this research. The objective of mapping relative permeabilities in 
the interior region of the three-phase ternary diagram is not achievable through 
dynamic water floods. The flooding of the glass bead pack at residual water satura- 
tion and several combinations of decane and benzyl alcohol saturation, resulted in 
the development of an oil bank and subsequent bypassing of the interior region of 
the ternary diagram during the flood. As can be seen on the dynamic flood 
recovery curves of runs 4 and 5 (Figures 7.51 and 7.56), the water flood displaced 
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the decane to its dynamic residual saturation at breakthrough. The Buckley- 
Leverett oil saturation jump at breakthrough, caused by the development of the oil 
bank, can be seen on the saturation trajectories of runs 4 and 5 (Figures 7.52 and 
7.57). Figures 7.52 and 7.57 clearly show that the interior region of the ternary 
diagram is bypassed. 

During the water flood and the development of the oil bank prior to break- 
through, two saturation shocks did not develop as may have been expected. The 
water saturation behind the oil front is continuous, in other words, at breakthrough 
both oil and water are produced. The results did not show 100% oil production fol- 
lowed by the appearance of water. Figures 8.2 through 8.4 show the core saturation 
profiles at breakthrough for runs 4, 5 and 8 respectively. The data for Figures 8.2, 
8.3 and 8.4 are contained in Tables A.14, A.16, and A.18, respectively. It may 
appear in Figure 8.2 that a water shock developed. This is not the case, however, as 
the recovery data in Table A.13 shows that a small amount of water is produced at 
breakthrough. The resolution of the graph is such that one cannot see the water 
Saturation increase until the dimensionless position of 0.72. The increased initial 
saturation of the benzyl alcohol may be the cause of the apparent lack of initial 
water production in run 4. Figures 8.2 through 8.4 show that a significant oil bank 
is developed, even in the case of Figure 8.4 where the water flood was initiated at a 


low benzyl alcohol saturation. 


The development of the oil bank, and the resulting increase in oil saturation at 
breakthrough in runs 4 and 5, enabled the relative permeabilities to be measured 
along a continuous saturation trajectory. This is in contrast to the water flood 
experiment at residual decane saturation conducted in run 9. This experiment had a 


jump in the saturation trajectory along the residual decane saturation line resulting 
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in the inability to dynamically measure the relative permeabilities in that region. 
These differences can be seen in the saturation trajectories shown in Figures 7.27, 


7.52 апа 7.57. 


Of interest in the comparison of the dynamic water floods of runs 4 and 5 with 
the water flood of the core at residual decane saturation in run 9, is the difference 
between the end decane saturations. The residual decane saturation arrived at by 
steady-state flooding in run 9 is lower than the dynamic residual decane saturation 
achieved by water flooding in runs 4 and 5. This is seen in Figures 7.52 and 7.57 
when compared to Figure 7.27. It is unclear if this is a phenomenological 


occurrence, or if it is due to experimental error. 


There are minimal hysteresis effects in the comparison of the relative per- 
meabilities at residual water saturation in the drainage and imbibition series of 
steady-state experiments and the drainage dynamic decane flood. As seen in Fig- 
ure 7.48, the oil relative permeabilities are not noticeably different in any of the 
runs, whether determined by steady-state or dynamic experiments. The decane rela- 
tive permeability curve differed slightly in slope during the drainage and imbibition 
series of steady state experiments as seen in Figure 7.41. When compared with the 
dynamic flood, the decane relative permeability curves maintained their relative 
shapes and slopes, but are positioned slightly above the dynamic flood generated 


curve. This can be seen in Figure 7.47. 


Run 9 showed some shifting of the water and benzyl alcohol relative permea- 
bility curves when decane is present at residual saturation. As shown in Figure 
7.33, the relative permeability curves of both the water and benzyl alcohol in the 
presence of residual decane (B and D) dropped below the curves generated when 


just the water and benzyl alcohol are present (A and C). This lower curve shifting 
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in the presence of decane is expected. The oil curve showed a more pronounced 
shift than did the water curve in the presence of residual decane. The reason for 


the increased shift in the oil curve over the water curve is unknown. 


8.4 Questions Raised From the Research 


Several major questions arose from the project's findings. Future work needs 
to be done on determining the reason for the development of the oil bank and why 
there is a change in the residual decane saturation between dynamic and steady-state 
water floods. Work needs to be done to see if there is some fluid flow mechanism 


that is responsible for these observations. 


The interior of the ternary diagram could not be mapped with dynamic water 
floods. A topic of future work will be to determine how to map this interior region. 


Steady-state experiments may be required. 


8.5 Implications 


The extension of the two-phase theory to three phases is validated by the exact 
graphical matching of the core-end saturation change between the displaced decane 
and the oil bank in the dynamic water floods of runs 4 and 5. Examination of Fig- 
ure 7.51, the expanded recovery curve of run 4's dynamic water flood, is an exam- 
ple of this validation. The y-intercept of the tangent to the oil recovery curve at 
breakthrough is —0.36, this represents a saturation increase of 36% in the benzyl 
alcohol. The extension of the flat decane recovery curve intercepts the y-axis at 
0.36, this indicates a decrease in decane saturation of 36%. The saturation changes 
are equal, the decane is replaced by the benzyl alcohol. In addition, the independent 
flow-rate portion of the theory is validated by the unintentional flow-rate change 


seen in the first two dynamic water floods of run 9. As seen in Figure 7.19, the 
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normalized pressure profile is not affected by flow-rate changes if they are incor- 
porated into the calculation of P. The inclusion of the flow rate in the normalized 


pressure drop, ®, is permitted only in the cases were gravity effects are neglected. 


This research provides directions for future mapping of the three-phase 
diagram. If the phenomenon of a self-sharpening oil bank is true, it will have an 
impact on oil recovery. It may be beneficial to increase the gas saturation in the 


reservoir prior to a water flood in order to facilitate the development of an oil bank. 





156 
Chapter 9 


CONCLUSIONS AND RECOMMENDATIONS 


An experimental investigation of the simultaneous motion of three phases 
within porous media is conducted with the purpose of quantifying and qualifying 
immiscible three-phase relative permeabilities. This research is exploratory in 
nature and intended to indicate trends, provide a background for further research, 
and generate generic data for use in any field or discipline dealing with the flow of 


three phases through a porous medium. 


It is found that the interior region of the three-phase ternary diagram cannot be 
mapped with dynamic water floods when starting at residual water saturation. It is 
discovered during water floods originating along the residual water iso-saturation 
line, that an oil bank develops in front of the advancing water front and drives the 
decane out of the core with piston-like dislacement, bypassing the interior region of 
the diagram. This will have an impact on the future directions taken in mapping 


the three-phase diagram and may have ramifications in the area of oil recovery. 


Grader and O’Meara’s extension of the two-phase Welge and JBN dynamic 
displacement theories to three-phases appears valid. This is shown by the equivalent 
Saturation changes of the decane and benzyl alcohol graphically determined in the 


dynamic water flood recovery curves of runs 4 and 5. 


Directions indicated for future investigation by this research include: 


e Perform the experiments, and determine the saturations using the CAT scanner. 
The independent saturation determination will provide the ultimate validation 


of the extension of the Welge/JBN theory. 





[57] 


Attempt to map the interior region of the three-phase diagram by originating 
the dynamic water flood from a point in the the three-phase diagram where 
mobile water exists. Conduct dynamic oil floods in an attempt to enter the 
interior of the diagram. Utilize steady-state experiments to reach saturations in 
the interior of the diagram unreachable by dynamic floods. 

Conduct more research in the area of the hysteresis experienced during water 
floods conducted at and above steady-state residual decane saturation. 

Examine the effects of a consolidated core using the same fluids used in the 
glass bead pack experiments. 

Introduce an actual gas phase, such as nitrogen, into the idealized glass bead 
packed core and determine the effects of the differing interfacial tensions of 


the fluids on relative permeability-saturation relationships. 
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Table A.1 
Fractional Collection Analysis Summary; 
Run 2 
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Table A.2 
Relative Permeability Analysis Summary; 
Run 2 








163 


Table A.3 
Fractional Collection Analysis Summary; 
Run 3 
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Table A.4 
Relative Permeability Analysis Summary; 
Run 3 
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Table A.5 
Relative Permeability Analysis Summary; 
Run 9 - First Dynamic Water Flood 
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Table A.6 
Fractional Collection Analysis Summary 
Run 9 - First Dynamic Water Flood 
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Table A.7 
Relative Permeability Analysis Summary; 
Run 9 - Second Dynamic Water Flood 
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Table A.8 
Fractional Collection Analysis Summary 
Run 9 - Second Dynamic Water Flood 


ороророоророоо 


DIS SIS O IS SOS SO SI O O cee dee о = 
GNO CIO m5 


0. 
0. 
0; 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 








169 


Table A.9 
Relative Permeability Analysis Summary; 
Run 9 - Third Dynamic Water Flood 
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Table A.10 
Fractional Collection Analysis Summary 
Run 9 - Third Dynamic Water Flood 
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Table A.11 
Relative Permeability Analysis Summary; 
Run 8 - Dynamic Decane Flood 
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ТаЫе А.12 
Fractional Collection Analysis Summary 
Run 8 - Dynamic Decane Flood 
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Table A.13 
Fractional Collection Analysis Summary; 
Run 4 
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Table A.14 
Relative Permeability Analysis Summary; 
Run 4 
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Table A.15 
Fractional Collection Analysis Summary; 
Run 5 
tube tube cum. cum. total oil water | decane | PVI 
oil decane oil decane fluid rec. rec. rec. 
vol. vol. vol. vol. vol. (PVD | (PVI (РУТ) 
(сс) (сс) (сс) (сс) (сс) 





O0 -] ON tA 4 шо t2 — 
слее ле ее е ос 
Seen 
бле очсо осо е 


> 
— o 
о е 
о 


0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 








176 


Table A.16 
Relative Permeability Analysis Summary; 
Ви 5 
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ТаЫе А.17 
Fractional Collection Analysis Summary 
Run 8 - Dynamic Water Flood 
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Table A.18 
Relative Permeability Analysis Summary; 
Run 8 - Dynamic Water Flood 
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Table А.19 
Run Core Characteristics 
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COMPUTER PROGRAMS 
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B.1 READPT.BAS 


10 REM *** CONTINOUS PRESSURE SAMPLE PROGRAM *** 


20 CLS 

30 ON ERROR GOTO 500 

40 OPEN"COM1:1200,N,8,1,CS,DS,LF"AS#1 

50 PRINT#1, "*0100VR" 

60 CLOSE#1 

80 INPUT "FILENAME FOR RESULTS";FILE$ 

90 INPUT "ENTER CORE NUMBER:";TYPE$ 

100 INPUT "WATER RATE IN [CC/MIN]:";RATEI 

110 INPUT "BENZYL ALCOHOL RATE IN [CC/MIN]:";RATE2 

120 INPUT "DECANE RATE IN [CC/MIN]:";RATE3 

130 OPEN FILE$ FOR APPEND AS #2 

140 WRITE #2,"DATA FILE: " ,FILE$ 

150 WRITE #2,"CORE NUMBER: ",TYPE$ 

160 WRITE #2,"H20 RATE: ",RATE1," BA RATE: ",RATE2," C10 RATE: ",RATE3 
170 PRINT "" 

180 PRINT "PROGRAM IS STOPPED BY PRESSING 'S' (CAPITAL - S)" 
190 PRINT "" 

200 PRINT "PRESS ANY KEY TO START 

210 X$ = МКЕУ5: IF X$ ="" THEN 210 


220 REM *** READ INITIAL TEMPERATURE IN DEGREES CENTIGRADE *** 


230 OPEN"COM1:1200,N,8,1,CS,DS,LF"AS#1 

240 PRINT#1,"*0100Q3" 

250 INPUT #1,A$ 

260 A$ = MID$(A$,6) 

270 WRITE #2,"INITIAL TEMPERATURE (C): ",A$ 
280 PRINT;"INITIAL TEMPERATURE (C): "A$ 
290 PRINT#1,"*0100P4" 


300 REM *** INITIALIZE TIMER *** 


310 NOW = TIMER 
320 DT = TIMER - NOW 


330 REM *** READ PRESSURE VALUES IN PSI *** 
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340 INPUT#1,A$ 

350 A$=MID$(A$,6) 

360 WRITE #2,DT, VAL(A$) 

370 PRINT;DT, A$ 

380 AA$ = INKEY$: IF AA$ = "S" THEN 400 
390 GOTO 320 

400 PRINT"PROGRAM STOPPED" 

410 REM *** READ FINAL TEMPERATURE *** 
420 PRINT#1,"*0100Q3" 

430 INPUT#1,A$ 

440 A$ = MID$(A$,6) 

450 WRITE 42,"FINAL TEMPERATURE (C): ", A$ 
460 PRINT;"FINAL ТЕМРЕКАТОКЕ:" ,А5 

470 CLOSE #1 

480 CLOSE #2 

490 END 

500 PRINT "ERROR NR";ERR; "ON LINE";ERL 
510 RESUME 490 
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В.2 ПАТАЕКАС.Е 


C 
C 
C 


O o 


This program does fractional collection analysis calculations 
for water floods rate. 
There is a limit of 50 collection containers 


dimension o(50), w(50), d(50) 

dimension Vo(50), Vw(50), Vd(50), Vt(50) 
dimension Vto(50), Vtw(50), Vtd(50), Vtt(50) 
dimension pvi(50), Itube(50), Ro(50), Rd(50), Rw(50) 


implicit real*8 (a-h,o-z) 


open(unit=3,file=’recovdata’ ,status=’unknown’) 
rewind(unit=3) 


open(unit=4,file=’recovprint’ status=’unknown’) 
rewind(unit=4) 


open(unit=8,file=’echodata’,status=’unknown’) 
rewind(unit=8) 


Variables are as follows: 
о = level of the benzyl alcohol in the container 
w = level of the water in the container 
d = level of the decane in the container 


Vo = volume of benzyl alcohol in the container 
Vw = volume of water in the container 
Vd = volume of decane in the container 


Vto = total volume of benzyl alcohol collected 
Vtw = total volume of water collected 
Vtd = total volume of decane collected 


Ми = total combined volume of fluid collected 


pvi = pore volumes injected 
Itube = tube number 
Ro = oil recovered 





со с с с с С с С 
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Rd = decane recovered 
Rw = water recovered 
Data must be entered into the data file as follows: 
pore volume of core - (porevol) 
number of collection samples - (nsamp) 
total dead volume of flowing system - (AB) 
actual decane flow rate being pumped through system - (qd) 
actual oil flow rate being pumped through system - (qo) 


read (5,*) porevol 
read (5,*) nsamp 
read (5,*) AB 
read (5,*) qd 
read (5,*) qo 


Determine the fractional flow of the decane and oil 
fd = qd / (qd + qo) 
fo = qo / (qd + qo) 


Constd and consto are the dead volumes of the decane and oil 
constd = ab*fd 
consto = ab*fo 


Initialize the Arrays 
do 10 I =1,nsamp 


Itube(I) 2 0 
o(I) 2 0.0 
w(D = 0.0 
40) = 0.0 
Vo(D = 0.0 
Vw(I) 2 0.0 
Уда) = 0.0 
Vt(I) 2 0.0 
Vtw(I) = 0.0 
Vto(I) 2 0.0 
Vtd(I) 2 0.0 
Vtt(I) 2 0.0 
pvi(I) 2 0.0 
Ко(1 = 0.0 


Кай) = 0.0 
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Rw(I) = 0.0 
continue 


c Read in data 


20 


30 


do 20 I= 1,nsamp 


read (5,*) Itube(D, o(D, w(D, d(D 
continue 


Ics 

Vo(l) = o(l) - consto 
Мо) = Уо0) 

Ко(1) < Ую() / рогеуо! 
Vw(D 2 w(I) - о(1) 

Vtw(I) 2 Vw(I) 

Vd(D = а@ - w(I) - constd 
Vtd(I) = Vd(I) 

Rd(I) = Vtd(I) / porevol 
Vtt(I) 2 Vto(I) * Vtw(D) + Маи) 
pvi(I) 2 Vtt(I) / porevol 
Rw(l) = Vtw(I) / porevol 


do 30 I = 2,nsamp 


Vo(]) + о(1) 

Мо() = Vo( + Vto(1-1) 
Ко(1) 2 Vto(I) / porevol 
Vw(I) 2 w(I) - ОСТ) 
Vtw(I) » Vw(T) * vtw(I-1) 
Vd(I) 2 d(D - w(I) 

Vtd(T) = Vdd) + Vtd(1-1) 
Rd I) = Vtd(I) / porevol 
Vtt(T) 2 Vto(I) * Vtw(T) + Vtd(D 
pvi(I) - Vtt(I) / porevol 
Rw(l) = Vtw(l) / porevol 
contInue 

пазатр = пзатр + 1 


write (3,100) ndsamp 
cl = 0.0 





40 


50 


60 


80 


100 
101 
110 
120 
124 
125 


FH 
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c2 = 0.0 
write (3,101) cl, c2 
do 40 j = 1,nsamp 
write (3, 101) pviG), Кой) 
continue 


write (3,110) ndsamp 
write (3,101) cl, c2 
do 50 j = 1,nsamp 
write (3, 101) pvi(j), Rd(j) 
continue 


write (3,120) ndsamp 
write (3,101) cl, c2 
do 60 j = 1,nsamp 
write (3, 101) pvi(j), Rw(j) 
continue 


write(3,124) 0 
write(4,125) 
write(8,126) 
write(4,150) 
write(8,150) 
write(4,130) 
write(8,131) 
do 80 I = 1,nsamp 
write(4,140) Itube(D, Vo(D, Vw(D, Vd(I), 
Vto(I), Vtw(I) Vtd(D, Vtt(D, Ro(D, Rw(I), 
Каџ), рукђ 
write(8,141) Itube(D, o(D, w(D, d(I) 
continue 


Stop 
format (I5,' Ro data’) 
format (f12.5,5x,f12.5) 
format (I5,' Rd data") 
format (I5,' Rw data’) 
format (I5) 
format (‘fractional collection analysis’) 





126 format (fractional collection data ’) 
130 format (# Vo Vw Vd Мо Ум Ма Vtt Ro 
$ Rw Rd РУГ) 
131 format ( # O W а ”) 
140 format (I2,1x,f5.2,1x,6(f6.2,1x),4(f5.2)) 
141 format (13,1x,3(f7.3,1x)) 
150 format (’ ’) 
end 


187 





188 


B.3 PVST.F 


This program takes the raw pressure vs time data and averages each 
"n" number of up to " Istop " number of raw data points for 

fixed injection flow rates. 

The pore volume and flow rate must be provided under the parameters 


с © 0» 


implicit real*8 (a-h,o-z) 


dimension time(500), press(500) 
dimension pvi(500), phi(500) 

dimension avetime(500), avepress(500) 
parameter ( п = 100, Istop = 10100) 
parameter ( q = 2.2175, porevol = 134.01) 


open(unit=3,file=’ pvivsq’ ,status=’unknown’) 
rewind(unit=3) 


open(unit=4,file=’ pvst’ status=’unknown’) 
rewind(unit=3) 
Icount = 0 
npoint = 0 
read(5,*) time(1), press(1) 
1 sumtime = 0.00 
sumpress = 0.00 
ncount = 0 
do 10 I = 1,n 
Icount = Icount + 1 
ncount = ncount + 1 
read(5,*) time(l), press(I) 
Sumtime = sumtime + time(I) 
Sumpress = sumpress + press(l) 
if (Icount .GE. Istop) go to 20 
10 continue 
20 npoint + прош + 1 
avetime(npoint) = sumtime / (ncount * 60.0) 
avepress(npoint) = sumpress / ncount 
phi(npoint) = avepress(npoint) / q 
pvi(npoint) = (avetime(npoint) * q) / (porevol) 





30 


40 
100 
101 


if (Icount .GE. Istop) go to 30 
go to 1 
write(3,101) npoint 
write(4,101) npoint 
do 40 Iz 1,npoint 
write(3,100) pvi(I) , phi(I) 
write(4,100) avetime(I), avepress(I) 
continue 
format(f16.8,3x,f16.8) 
format(16) 
end 
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B.4 VARIQ.F 
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о Сс су с; 


о 0O 0 CO 
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This program determines the normalized pressure profile vs pore 
volumes injected for varying flow rates. 


The program first does fractional collection analysis 

calculations for water floods then calls the subroutine fixflow 
which determines the normalized pressure profile for varying flow 
rates by determining the flow rate for each collection time. 


There is a limit of 50 collection containers 


dimension 0(50), w(50), d(50) 

dimension Vo(50), Vw(50), Vd(50), Vt(50) 
dimension Vto(50), Vtw(50), Vtd(50), Vtt(50) 
dimension pvi(50), Itube(50), Ro(50), Rd(50), Rw(50) 
common/a/ ctime(50), q(50) 


implicit real*8 (a-h,o-z) 


open(unit=7 ,file=’datafrac’ ,status=’old’ ,access=’unknown’) 
rewind(unit=7) 


open(unit=3,file="recovdata”,status="unknown”) 
rewind(unit=3) 


open(unit=4,file="recovprint”,status="unknown”) 
rewind(unit=4) 


open(unit=8,file="echodata”, status="unknown”) 
rewind(unit=8) 


Variables are as follows: 
o = level of the benzyl alcohol in the container 
w = level of the water in the container 
d = level of the decane in the container 


Vo = volume of benzyl alcohol in the container 
Vw = volume of water in the container 
Vd = volume of decane in the container 





e 


ооо о 


сб со с) сб С с) 
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Vto = total volume of benzyl alcohol collected 
Vtw = total volume of water collected 
Vtd = total volume of decane collected 


Vtt = total combined volume of fluid collected 


pvi = pore volumes injected 
Itube = tube number 

Ro = oil recovered 

Rd = decane recovered 

Rw = water recoverd 


Data must be entered into the data file as follows: 


pore volume of core - (porevol) 

number of collection samples - (nsamp) 

total dead volume of flowing system - (AB) 

actual decane flow rate being pumped through system - (qd) 
actual oil flow rate being pumped though system-- (qo) 


read (7,*) porevol 
read (7,*) nsamp 
read (7,*) AB 
read (7,*) qd 
read (7,*) qo 


Determine the fractional flow of the decane and oil 


fd = qd / (qd + go) 
{о = до / (qd + qo) 


Constd and consto are the dead volumes of the decane and oil 


constd = ab*fd 
consto = ab*fo 


do 10 [= 1,пзатр 
Itube(I) 2 0 
o(D = 0.0 
w(I) = 0.0 
d(I) 2 0.0 
Уо(Г) = 0.0 
Vw(I) 2 0.0 
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Vd(I) = 0.0 
УКГ) = 0.0 
Vtw(D = 0.0 
Vto(I) 2 0.0 
Vtd(I) = 0.0 
Vtt(I) 2 0.0 
pvi(I) = 0.0 
Ко) = 0.0 
Rd(I) = 0.0 
Rw(I) = 0.0 
УКГ = 0.0 
q(I) = 0.0 
continue 


do 20 [= 1,пзатр 


read (7,*) Itube(D, o(D, w(D, d(D, ctime(I) 
continue 


EN 

Уо(Г) = о(Т - соп$ю 

Моа + о(1) 

Vto(D = Vo(D 

Ro(I) = Vto(I) / porevol 

Vw(D = “(0 - од) 

Ума = (0) - об) 

Vtw(I) » Vw(I) 

Vd(I) = d(I) - w(I) - constd 
Vda = d(I) - w(I) 

Vtd(D) 2 Vd(I) 

Rd(I) - Vtd(I) / porevol 

Vtt(T) 2 Vto(I) - Vtw(I) * Vtd(I) 
pvi(I) = Vtt(I) / porevol 

Rw(I) 2 Vtw(I) / porevol 

Vt(I) = Vo(I) * Vw(I) «* Vd(I) 
q(D) = (Voa + Vwa + Vda) / ctime(I) 


do 30 I = 2,nsamp 


Мо( ) = об 
УюШ + Мо) + Ую(-1) 
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50 


60 
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Ro(I) = Vto(1) / porevol 
Vw(I) 2 w(I) - од) 

Viw(D = Vw(D) + vtw(1-1) 
Май) = 40) - мр 

Vtd(I) 2 Vd(T) * Vtd(1-1) 
Rd(I) 2 Vtd(I) / porevol 

Vtt(I) 2 Vto(I)  Vtw(I) + Vtd(I) 
pvi(I) 2 Vtt(I) / porevol 

Rw(]) = Vtw(l) / porevol 
ҮШІ) = Vol) + Уу(1) + Vdd) 
q(I) = Vt(I) / ctime(I) 


continue 
ndsamp = nsamp + | 


write (3,100) ndsamp 
cuu 
с2- 0) 
write (3,101) cl, c2 
do 40 j = 1,nsamp 

write (3, 101) руцј), Ro(j) 
continue 


write (3,110) ndsamp 
write (3,101) cl, c2 
do 50 j = 1,nsamp 
write (3, 101) pvi(j), Rd(j) 
continue 


write (3,120) ndsamp 
write (3,101) cl, c2 
до 60 j = 1,nsamp 
write (3, 101) pvi(j), Rw(j) 
continue 


write(3,124) 0 
write(4,125) 
write(8,126) 


umm 


"A > 
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write(4,150) 

write(8,150) 

write(4,130) 

write(8,131) 

do 80 I = 1,nsamp 

write(4,140) Itube(D, Vo(I), Vw(D, Vd(L), 
Vto(D, Vtw(I), Vtd(D, Vtt(D, Ro(D, Rw(l), 
$ Rd(D, pvi(l) 
write(8,141) Itube(D, o(D, ма), 80) 
80 continue 


fA 


call fixflow 

stop 
100 format (I5, Ro data’) 
101 format (f12.5,5x,f12.5) 


110 format (I5,' Rd data") 

120 format (I5,' Rw data’) 

124 format (I5) 

125 format ("fractional collection analysis") 

126 format (fractional collection data °) 

130 format( # Vo Vw Vd Мо Ум Уа Vtt. Ro 
$ Rw Rd PVI) 

131 format ( # O W d >’) 

140 format (I2,1x,f5.2,1x,6(f6.2,1x),4(f5.2)) 

141 format (13,1x,3(f7.3,1x)) 

150 format ( ’) 


end 
CFE Rk k k ak ak ak ak ak ak ak ak ЖЖ ЖЖ ЖЖ ЖЖ ЖЖ ЖЖ ЖЖ ЖЖ ЖЖ oe 2 A A ОК 


сжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжжж 


subroutine fixflow 


c This program takes the raw pressure vs time data and averages each 


c "n" number of up to " Istop " number of raw data points to give 
c a managable data file for the graphs 


c This program also adjusts the pressure time data to account for the 
c delay in starting the run and the dead volume of the core 
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implicit real*8 (a-h,o-z) 


dimension phi(2000), pvi(2000) 
common/a/ ctime(50), q(50) 


set the pore volume of the core: "porevol" 


set the constant displacing fluid flow rate: "qf" 


set the delta time (seconds) from the start of the data file 
until the beginning of the run: "deltat" 


set the flowing dead volume of the core: "Vdead" 


parameter ( porevol = 131.97) 
parameter ( qf= 2.2175) 

parameter ( deltat = 103.0) 
parameter ( Vdead = 1.82) 
parameter ( п = 100, Іѕ‹ор = 19400) 


data.in contains the pressure time data from the run 


open(unit=1 1,file=’data.in’ ,status=’old’,access=’ unknown’ ) 
rewind(unit=11) 


fig.data is the output file that contains the corrected normalized 
pressure and pore volume injected data 


open(unit=12,file=’ fig.data’,status=’ unknown’) 
rewind(unit=1 2) 


calculate the time adjustment: гад) 

tadj = deltat + 60.0 * (Vdead / q(1)) 
set initial values 
Icount = number of actual data points read 


Icount = 0 
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Npoint = number of points recoreded for graphing 
Npoint = 0 
ntube = collection tube that data are for 
ntube = 1 
runtime = cumulative time of the run 
runtime = 0.0 
tubetime = cumulative run time that the tube "ntube" was collected 
coltime = ctime(ntube) 


Vinj = cumulative fluid volume injected into the core 
corrected for Vdead and time 


Vinj = 0.0 


Vinjole = cumulative fluid volume injected into the core 
at the beginning of the averaging period 


Vinjold = 0.0 
Vave = average volume injected into core during a time period 


Vave = 0.0 

rtold = runtime at the beginning of the time step 
rtold = 0.0 

qave = average flowrat during the time step 
qave = 0.0 


ctime(I) = collection time of tube "I" 

q(I) = flowrate of tube "I" 

told = last time point read 

dT = time difference between last time read and current time read 





999 


1 


C 
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do 999 j = 1,38 
continue 


read(11,*) t, p 


told = t 

Icount = Icount + 1 

t = t - tadj 

if (Icount .GE. Istop) go to 20 


if (t .It. 0.00000) go to 1 
Npoint = Npoint + 1 


phi(Npoint) = p / (q(ntube)) 
pvi(Npoint) = 0.0 


start loop to determine pvi vs normalized pressure 


sumpress = 0.00 
rtold = runtime 
Vinjold = Vinj 


Vave = 0.0 
ncount = 0 
до 10 I2 l,n 


Icount = Icount + 1 
ncount = ncount + 1 
read(11,*) t, p 

dT = dabs((t - told)/60.0) 
runtime = runtime + dT 
told = t 


if (runtime .gt. coltime) then 
ntube = ntube + 1 

coltime = coltime + ctime(ntube) 
goto7 

end if 


Vinj = Vinj + dT * q(ntube) 
Vave = (Vave + dT * q(ntube)) / float(I) 
qave = (Vinj - Vinjold) / (dabs(runtime - rtold)) 
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sumpress = sumpress + p 
if (Icount .GE. Istop) go to 20 
10 continue 
20 Vave = Vave +Vinjold 
avepress = sumpress / ncount 
Npoint = Npoint + 1 


phi(Npoint) = avepress / qave 
pvi(Npoint) = (Vave) / (porevol) 
if (Icount .LT. Istop) go to 5 


C print data 


write(12,100) Npoint 
do 30 I = 1,Npoint 
write(12,200) pvi(I), phi(I) 
30 continue 
100  format(I5) 
200  format(f16.8,3x,f16.8) 
retum 
end 











| WERT 
| BOOKBINDING $ 


Qrantville Pa 
OV-DEC 1990 


N 


e 


с дна Bowne 





Thesis 
T4174 
el 





Thomas 

Three-phase dynamic 
displacement measurements 
of relative permeability 
in porous media using 
three immiscible liquids. 





















































































































































































































































































































































































































































































































































































































































































































































































































Ена 
РО 
dO rr 
Me РАР 
m ДАМА 
т ИЕ 
x Tu Ж.А. Bore er vars 
PES tate 
хе (и Y NI ната 
c T "gw me 
Се P i 
1 АНХ vá ow SEN 
TAN е Наро ма 
Е И a Om: + #6; A 
2. NL 
| игра Ms Pini 
n Hin К br re Хх н Jn “ [1 а? 1g 
ut | ee Е 
Дои сорта ЖЕСЕ ида КТЖ ЗАН 
45454 М ЖЫЛДА fae я о ju è T Ad e gs 
f AE „он а Дон M mas IS ДҮР пря E (БРАН 
‚сй ita 153 ей и ФА se ни wi ud EN 2. 
TM 2. Е DEN үй ДЯ пен [4 
E с T n, E м pn СК D PED PM 5 
И By Hp ес 5 it A. v НЕН 
Lee pur feet m Ai ae ЖҰ КАН 
СА I Set PL M AI ду Liga UT И Ba ұя 
3 4 А В Б а 5 RR таа м te pti "x "ura 
АСА ДЕ BER H КЕН er TAE ios n d s dona На 
Є uS "P д Tm И COUR Жа Party "uu Py En не КШ f. 
m де ара ат : rere En о P ULL TRUM Pet a Mi rn "i. k 
е; 1 К mi =. Е аў, к | a Е $ МЫ ы ААА ДНА 
Ун ie QU си Зи Бара Е. N duod 
m эл Се Hd 7 m FE ге "3 E más iR Е, pb И: Ди bra ы D ег А TY (5 Dl $ Ки 4 
m. des а Ден er г t AE АУК 4, А: Ве n ny LU s МА TREO 
xo де и к — ма AC ups 212 
* = > H am ДА Ohi 27 ieu Й : mon ПЕЧ fecal 
ria m Е ie pis и D Ме CANTE м RO n d AS Ни м PW NE: ТК 
T Ва АУ e Ba oe jen . 4 
ШУНУ А. dar ДА y “ RERIT Ame уђе У На НЕ жу Pl Ц TN P m 
ru UT шен суса: vut 2. HM А. at e И Doe RE ІШ 
m A сы м Fede НА att 242 t den як a T . m TH UL ¿Ls peri уда ЕД: xu | 
irá P 1 x» НА E b 0 т. КИ Ра 3 
Hase КОИ; ee o | НИ 
дет rep ўа. 244% i eh И ЖЯ КРАН | ин + | 
EU, ой 2 ЊЕ ! ||| 
д; г -Ф 7411 ки $ { К ЕЛУП. яб аы LI DAN H i AA n И jur “ |". #4 АМ | | ||| 
552% Е TÉ ИЯ а r Ыы He nam ті; uu xh QUT. Наз IET Ц УҒА 4 м Нео r t НФ ) | | 
m Dr mc Я EDU. К е?” ЖАҚАУ Aa il ie lage Й ТАП n PN" JM UT sr. 4 nU E ДА " Ги Не ІШ | ІТ 
T naeh. Ван Ру: Да Мег At 26 JA пи er AS hi in 7 н ЕА ДАЈЕ те нн ||| 
бити ү с Сори er ER nen m пр n no аз, nn | | 
ut ЖУ ХИ p E а Ш ИИ 
ө (4 Lr d ДА H Р a ? ana и m ДА n ur КЕПП ЖЫ РИКА. Ж | | 
d 52.7 НА А Фр Ehe ГИ mr Tm Х T 51 m И) br «е: л D a “ Й ў fi eel А Ti ‘ "Nu 
; Р yr T PU ы” о кп РЈ ЦН Dr h: TI M AP | 
? 222 MT 7243 “6 ра e PM Fe 4 bo m 4 я ‚е ВЯ T а 44 gh n Tu E ШЕ, МА AE КЕЙ Pit H | 
Да М б ТІ p pe = Mae 2 иа 274 РАД RH ELIT] LI] she T "T A E EU EJ SACRE E Й Riy or Ц -t | || 
Б = H D Ser yd is Aa ГА M ped И) ОЕ ЫЫ 2, ' РЕН КАРЕ [| EN “ П [4 PAL [ i || 
та, ett NT: PE Mire мы n JI И ген TM Ha ши Е. г Ap ДА PERT nm Eo бет | | 
re и nje han И ME e ИД. Қы Хе par ees T v ped в parece ДА ИДЫ олы жи МАЈ ue nts | 
gi а arg; Rt АН Tn 2 =, Hr и Дед: ии Р; +: 4 ПРИ ня ЕЖ. В Mid Н. 
iets OSA Pd ns mu eS О пи э НИ z НА dnt Ран Я Кр ДИ к у рено Н 2 ҚЫ pu y и ? || | | ES 
аи rj А e 5 БЕСТІ dien jt ne aS а пин MEE ЦЕНИ 1 dM "n i т; у 
га вом Ё H nd he у 7 4. КИ ^T 2.4 A ж. РАЛУ, HI ГИ ГА) hs ж. Pe О te E Г СА Да TY " Hi "T! "n П it ГЕ А Tn " МИЛА E Е А . H > 
"$ ARE AR Am M в pos ~ СА КЕ y m Bon aid К , П РНР МА аб 4 ГА “. 16, hs m A ( m LU E 24 7 Figs set “ А T ' ”. Н 
а ар , EH 4 T4 ae e ер е Ше Sae i JM Aa A EY God Ent бы НА пра MI um PPS 2 NA iud га Sy , 7 АЎ vt ы 
E уз о BE Пе 22. по e DOS | AS (cM 
ew ELE VE LEE и | ы К ae te PA RID yu Ақ Ни D Џ ec AD A Ц "i A. n e е. DU Pa. OE nu P rU 
< :2 ed "E Fer ue А / aud 4 24 И n» ДАН и С и ИЕ HER 31 uu ЕА at ае PR А JL ВЫ u. АК Ж, re SEES жанын 
P Er ig I» > n РИ 596 Е 42 fe 2 "v В Я "T ‘ pod ! e. nues , t et d йе á te tot m беа ГЕ Й м Ку t t гы fi А: D “а 14 g " ? и П [1 , 
Ps СИЛ FE F nn PHA 4 ЕД Пее 2. = dues P E mu UU. КИ us Қ, На a o zt r E қ .: i я А nt Hr Е. А Ax "i ЈЕ р E LI 
Мера ard ada Al я AR РН Фара e б? Hd vr ne Мин они UL HÀ apu Ане et qn ии ДР. d AP ta s Пк Pon TM "D «ДАЛ ti Paa Ја Ж ix 
VM Я АКА ма ДАЛ МО ў T с [Л "C пе a LE т Ра Pr) Pod ТҰ ХАЙ " Жала NI 7 m 7 АЈА РИА 1 пре T P ( Pun 7 А fo ; в а Й П LI 
Е А T rå ДРЕ РАЗ БА "і И e I Nf er ж ГА 5 ПР й e... Ф. р " ЖІ) м 4" " E n ote E ¿1 8 в Age £ [D VR " D ш on " 2 [E | T И П (a ' ' И р qr 
92 E TEN m - + НИ БУНЫ "m ги 42 ё. й та беду n m (17 И) er” М ДЕ “4. .. г ШП E fe T mi i ШЕ а PME t r " Н "E » ee te u , " ! E " D L П Й ! .. ГЕ LI " Й a 
2 ғ d ae E P прай ЫНЫ та Ы FAT, Жи d M РА NS ME Я УНАН ir gre Fad mU Un Pr. n] е Му ет А ' Pom nis На ы М ns На. EIS СЕГІЗ 
р uL p Е RER А р Е > крои са пе пий на er М VEM TOM 
ж” D E JF ES T" RR к i E "ET | ие RR] ee қ ЕЈ р bé as LE Д de + Мы ЫЫ М a A Ave as chai Й Ц mo ve єз М uL + ы ЫШ А. n E " REL П Да р Ц it AM ко. * Қ. Ы ] 4 eg 
ред © УИ 6m 205 ии РУХ ыр pi fd Pat) an Ei а ЗА ЈАЈА L АЕ ЕК НАД ЫДЫ еы de ‘ и а.” Ди ма TN m MN у ЖЕ FE у 
А у о E ова дау, ^h m A 2221 2 n i Ра u. "M Ex и se 3 ы "И ARE oS i И J ; a bi E е T ы с и 50 ки 
P ЕР D 24 ү a "T. "ё. 26 „ = Е H Р s ГЕ er} Ы A i a ААД Д Ф ИП ee Н Й 2 t Е E " ш " H В ЫШ 
BER пара => РАКА о г. . c TM ro DM E D ptt оН oe i DO р ы . He н T v. = oa р а CA e m 
ФАР, ма US Е 4 ~ Pes Ée ди АЛАНИ АНЫ UE UP h t EE y "PT АС > m ro: Ы NAA h ба ср: TES 2 ngu P d Да T Boy АЛЕ re е) ді. "i D що) ex аа 5 ра RA 
35 [943 n P n ff ej fide ЖУ ВИ КАП m 2 HT é E y п HUP tee " ’ ” N A йо А " м m ! $ 4. ¢ " "me " , 4 П .4 П О и de 24: П PT РАЈЕ La “ Es Й à П Ц E П m П | D 
33 E ty, е ТЫ s EPI H ВЕТ а П Ц nn Й Й РЕ “ КЕ E Tu ot P O " р h. М ш AT ee П Це $ '. [15 В Я t mE r ДАН M P y > E $ ! П у E E 
i E ry DE FEY dite dir АДИ ҒИ ҰМЫ an SNP ONES: EM us 0 NN pre Ы EET о Й г 517 ы . f 
UE: rr F : РА таз QUIT e , E E А МИ 4“ ДА "Pe ВЈ E d 0 E Ti + 1A E ‘ of ER ' ' Й П 2 a Ц Ц p LED П 1 AM 4 5 D Џ D P r . ы 
Y p dE РА НАДН BEER Ley е Над an MEM 17" pel EAE ПО и М A Џ 4! D p f E a € у , р + ын * E П FS Y E ре E e и В р 
44334 ПА n as Р) Es Ft й ы n Е Po Бы AE « г; Ш t Pee VM ES TIS ы L a ы ' Ballen г ' d Ы 
n BEE p er А Вие да TONO TR) 2 E i s RET Nd Ru do S pos EU RE А pants 0 n "e n 
ЕТЕ is, ХЕИ Ir Ед N ' mor zu. Far Br ge 5 “ г Р T nee П pn^ "n ЖЕТЕ u. ut ^ h m 4% Hi ^ е Ма 4: "n А, А г P 42) я [ "E E MEE ы я i re АЈА а Й y ў р 
ia d IUE ás бу Е я Tra M Ten mA б Нуи А (ЖА on gu ' E НА we ee A tai x jt PE "y En EIN e үну PM "P d e K ) E E SE иа wa 
t. s oy fa „82 » “4 62. A Wie АЙ AE E НАЕ 2% Му КАЛИ P dyes é Lou I, n E и py d M р TE 3 A P ^. Д " DL Rn қ Я Ba 5 D а И E P D 4 А = t И, A и P gi П Ка: Ри аа сш; : | E T А uw 
Hi АТ У FO "DA RE ЧЕ ШЫ; Жақ E py ия ДИ: AS DU e e и. 3 Ж Е КІЛІ А НЕНИ ў > Ж. % и D кх . т k 1, р ры ni и : * да ся = E D , М 4 д; Е, D E E , 
Ба, 2 с ба ИР E nos Е 5 ; : | 5 Ы Ш " . 8 р Ы 5 E ue ' M. М р 
EN A ta ee Usi LE it ЖА Гам ? Ж б, “у, D "i ND 2 E : e ПК Fd y и n Ру i R d eN n 0 БИГ: Қы, у 4 ^d "a e ' Ди р м м Нел | ae К у 
аи 22% E : PS A ri КІ ay КНИНУ пе M . и PL G м Пера ie es : DOCU d c d x y БАЖ, Ре ! A. 
Ея ХАНА u e einn LE EZ D HI Ы К 4 ғ ИЕ М ARM М et + Аси r Ц Fir er D Е H P Ы D e t Ph [ у у р E Џ 
3 dus 2 H P os X 47 gue bie А > .. 4 А Й oe " Д О .. г t Е 1 П E i " " “ " Ц mr . П Г ща. ot 4 Й А . Џ А " Й М " М р Й ] П 
MONA 3 б p en Ta jur 7 * Ц РАЯ EM П a; A n г a E D Й hg E Hn D m ЈЕ ты. ы А y р E ҒЫ у Ви : © E 
У ИД“ jr zd а ХО nd KS M и nan Ry ас a (ку р MUS Қа Да oon DX yos m T ot Lag DE А us i И a. 
Бра o H Ч И Y 2220 y i '' m и 5 г or y x н ЖУ) DOR р neu an у VM T i E Re A TEN uu * ша TEES р j 
3 АУЫ Ы ae Е ' А Й А i М у P G е Қ А 
де 3:4 88 As ү КЕ „ e "s m “ m Py bip Un ^ E са Р КІЛІ; 2 s M" = КЕ. о А. EE EL | Е | А : 
EA ЕА er Do 3» ou LE ter SA PR ede Rus “ "t na P au D D g 
МЕ Са Зы zh і ы ГЕТ П .% ыы . М и 9 шш PE ' Ы T T E ы z П ы кү Ц Ж ae Ы П ., a П * 0 E А М M . О t En LI П 
t Xe fu Й .. 0 А 74 2 С er , . aa " D А: Й П ‘ Й АЛЕ, # А Й те А 0 % П .. E " " қ "m M я П А P 
ПО и ” т. " ..: P rus си RE DO " ED Cx M N E a " П Ы О 
И D ҚАЛА | n D , en PEE : P. н кыр, ч 
Ц уа Чез NU Й DIM 2 DD e а, Fu 5 id ou у ary LUE n Еда n Ба "T caa o Џ t X .. 1 ы ЈА. и б x TRE E 
Н ЖИКА и LEE on DIU M: PUE DNE D EDS. estt, LL а ИН E . | ; rans кш 
DES Hd а Ји P H nf utes ты DD ы ү ш: "9 B | d ER % О or Ms " Ré ^y i. D = 200 Nx nh OU uo m A а t | | m ЈА 
Да Dre: MA Оо Пы n Й Po DO LEN D Ц n Й Й L -. , Џ Џ 
wy ys Je Най Да ODE A Й ы i DOE ы LX р а. Џ я d te ы Жы BEER: М ЕЈ ДАНИ м d e a П d 1 у d ' i р .. 
A 4 Е ш "PS ist 224 » ; ТИН N p ed ai E a Ха "EL S и М a Bd en ES SE о k ЕМ " m By "i ©з ч у P ТШ бъ: E ta П Ы Б A " ~ i Ш М , % у y М 
M г ч ВАК - УИ. Па ИА Ера "e и ТЕК А A oa tu ҚЫСЫ FE "us z S eN TON " E. AE А SN yo e .. 2 T аа "EE А пали „1: я 
зря Un та +, 7 зе j- ы "d КҮҢ КЕ ОР es. к Й е М 52 na О а СИ Ц балы р да РЕ D ME П А Жу qd Fark d E de D Hu e џ Р 4 D MU » D D 
иа Дер Щи КОШ ғ Дара HS Й 2t RS АДА = На ды DE На у ы А П E BA А ҚЫ ЊЕ ТЫ y а au a nas а "n Ey и ' "EL e қ "^u и, Й Ш . 
Ни 44 M E Uu 7 44 или na in: 7.3 E "a 5. Loos non ге b у Ва A ee a: при и Де b PEL s n d ^ NI ig A ки + A а T "n "i АШ m M у 2.4 
ri > г 2, "IS к - m ои rn р Ы РЫ E К ДЕДЫ "T | M р а P Re у у E р Р | қ КК 
M d ја о Do и A УР е DOO rr Ху ЗА 3 z В DO er us В ә КАЗ 1 a у 00 С T E б и Џ | Џ "ut у М ID 2 ve y Қ, x d x Џ 4 
dA gud г Hn LM Де ner ARUM Ы най, m "n ЖЫ "b "mu М P M ' E Ba | 
M е реа Ц ада у WA "y 2 DE Ex EAE d ғ не y ee T T и г AP PR Ы VU М е а i y = м un Й M Ar П А n. Eo M B Ш X ' D ү и Ы У n Ы Жы ; 
SES sd Lb Б, s е ul Ке Ы En ER К x pn и Ше Е "i Д 3 Ы p i 2d е Kan um s DA AN "n n s 5 e n " un j^ Ц и Џ Sn 
+: evened rou. е UN EON ei A Оа a TS А ве 1: Bars " О а О 
Е Y 58 rem DL ES ar е "n H РЕЛ Paes wd Mo ee о DD L] MEI aoe пи | П П A "m И а à fi M В бы DS oe И КК В EC Sg А 
443 ДЕК INE А ЫИ ee 2. mn "n Eo Касии E BE у қ у EI 
A y: + г Dee В 9976 ы , ы .. i PP .. А T ss па D .. 4 қалы; A es "ED L Ba of Ц x UR Ц .. D Й L PEN А у m ш 
БОЙ ] x :0 DU e Й =. Н Е d DN ' m DeL. qe d TET DOO Tod Ве ! Pane Hans a. D Y at ПЕ р ы "de 
dA UE 4 М . ч PS А “ 4 e ; ae .. FM ТҰ M DD P LL Е E Й E A © ДЫ $ E .  . " M И i Е 2 " "i " № Е гі; А А А, СИР ' d А П П t у "t б ' 
Stee ДА d Ped V MM D к E о Я о Бет Да и H ша 
Е 5 uf у Р ти ЛА Н E "LI 4 при 2a " i P , A DA i . LI "i 0 
0 A АА: Н М A f. T "m | зе BUT i ц% Pr d | EE x r 4 ys EN .. = ^d d Ни Ps , БЫ Ы Й + О "T П nmn зи ^u "T П of 5 Е Й i Ш + а L ҒЫ РЫ О | mn E a ' в a 
TERN Е "7s ir ЕА w" ғ Зб ik = Bit x. > .. $ Й аин dz Й .. 1 да "EON A ^ E .. А Жы И shire Қы P А Р " ws Е L D ^ 
rar Tg а Dr TS ar E АЫ à cu HE М 3, [] * NS у Ы М P па Las 
"te = " E E Po ~ 4 г E ef E и НД 4 DS n n je oe ии ic sd у Ў en Шш no B Е Р 2 е " 2 О .. и АЧ P А P О 4 Ц А Фа м Е ы А & s a 
Py E Tl г на 7 1 > = ^ ITE e. iy ee . E КТ м ? А ° Р . it D at А " toa А Б a 5 DE " rer Pr Й ad Е . ТЕ r А П . .. б f Р О РЫ И Й a 
Уфа d НЕ é FIR Ни A Adie Ld Й Ра > е. y Li 41 L Ы Ц [n E Ц ве LOO 4 по р О Ц en ы 0 Li p. р Ц : б, Ц h P қ ы Ы 
Зе 7 He 4: = А H ў У m 4» Ри H я " te М "ME ‚$ ot . г d ка е .. , ELA i i бе У RU d * P О 0 ы LU кын ‘ Ц * ы "d d 0 E б М с dos и U Ш D Ы Fig [] v 
4 ки H^ the: БОН 2 таи их, На . A Ex P қ к a ; D utt . NS D улу р Е E P "d ш oT has i 
14. y FI 5 ҒЫ 4 XP 1 E А, € E ВАН E ne n» б 2 RR за " ' n А oe P А ТЕ ГЕ. 2 n a И А "hu y P В РЕ К L ' Хы” МА N А Е de 
HI T2 t | MES 3 x , О Ha М net Й А да „* Ы y] ' О О . р d M А " 
y Ху. и : у ы 5% ҒАН ел ' | а ЕН ТА noo d TN Pe ae рад За ag ч ^ч x Ш x "i Жо о Қана: n D e Ж ЈЕ Es OS 
ч > А Ха D EXTA " ACE ғ У: , нощ ft Ц Т М р Abe 0 Т 0 " ' " 54 П Ы 
Lr DS ү еее ДАН . J^ ы Құр АСЫ У О Д z E NE ‘ i t қылы Е "P МИ y TE 2 i An n В p М 
БН ГЫ Й L4 L E r А i N Џ LN 472 О Й р "E D A D И ' 
ыы 3 ИЕ г 2re PIE a6 A EET мощи М LM р | ' iN у aon Зи р . 
È У! E ve Hd ~ М ГТ НКА Й бый BTE ES Er une б а oe AR ear а А m. р Pans 0 PES al E р ! Й И 
+ = LM (i FL А 167 e. А P А А г i " ща " у А P A А 4 x О А И a у " Л " e. n Ц М " 
1 е ? E r Á : LU. RC y ew С .. ei E Ч : у es ы ы О у 5. FUEL ш f 0 ы Pod D s Е, р G 
я y oe к Па “> КИП ПЕ: Rus: D ES Lot n А M ЖЫ A s 
"bi Ч АҢ В m yb Да е m Lm а. s ба тз s Ы ОД у они PN * да PX Td n . n и у я ub 
Eon =: Y ы ү ЈА, H are i T E U к "x t an Н x ч Ц Ву d ne 2. Да И а E E P А uL 1 ' i n "a 
YS XO : Jr р ve. a НАЦ ce = Ы ; Де MN ии ба ! ON 2 » mn Ц У | "s Ы me à i ма Ы D 
^ r ч Ур М E 8; £ ot 1 О за Ы FR Ш d 2 П , ? ~ ky E ы Е i 2 М ки f Ы oe 52 . НЕ ' Ы S bn d ч à т ы = "d M E Ру L a , "d М 
Хе: а E 7) В ¿E Dee пи е wo ГЫ Ы и, ' М Й я О ы Е Ц vo 
маза ПА 1 d А А Rod 8 Po A „+ ene ш oy Й А x tnn Пи " A P y к) 2 з 3 б As K L ҚЫ . А y O E "E А Ы 
де Би h in. M E А T» m А 5 А . i^ .. A P П .. E Ы c u“ 7 D Бы 9: « . Ц s D a " А .. Й О б, О > ы Е р 
~ y i Y БС ; E EE "E Ы 5 р G А Е о А X ы р 
Br М ЈА : 3 Ае Ха , ы | ES 3 i " s E P G Р р ubi Pi А > Т "d E Н Lo à ~ ; и 3 EN Е; в | Ы у d б. y у 
» На d н В г Д м Й р " D 
Patel derum e В t? Қат р ки. 6; ` ' ү DN DR ' пе "S ЈЕ ' in ee ine ae | à i 
ME Ка че 6 1 М КН ш E А М Ц П 2 Ч ҒЫ я ' ae и и Ы | у Ц Ы 0 Ы Е ' 
а Ч За n і Ы = О Ts ans? a. з L] E: ; я ... А е 2x а ê = Р Е À us Е > В 
низа. ~ wi ep QS у :\ Н К E ЕД. “” = и ^ . И пе А е XE Й тев. Р ОИ A E E in: "EM T 29% А 
a xo КК: мъ НЕ eae” Вн ЖАРЫ eut AD ии ғ Py е oues TEM y D г : 5 EN " Aa D ы ы 
Ir pre ; БЕ la ~ vt СИ ezo ta A Kr Ы A О A "m " . 5 Bi Е Ы SR z a 0 А, А се x " ча а Ы UM Ы Е 8 ВЫ но ES ПИ А 
1 ^m а ч +, 4. 2s Ц ae 0 = .. +. • L М Й 0 
А | тая базар сива Фа P. M ЕЖ КЕЛУ PR И a ЖЫ 2 Я, "on ҚТАР E M ed n EN n р t | a қР i "n pee р =. 
— at, Rate AN pur ыр Зи EN ne y . ат "s 5 e ey аи а E. i a Pian i ho ~ В р Е D ү 
H y er Б ы” n 2 <. rat ғат ЕА ща ; А Ms ы Р А ы we H Й 4 d of d Б E М Ц Ы 
етсек % — ys .. с Д c d 1522 E Ш M р d M d d 
» E A нү ВА 2 PEN = > ci " 5 0 = П П е А "E О : Й 
rS Le S Аа od: ‹ Ж x Em. E FE 5 P "d vat Ai Да T 2E. jJ "A А ES e Ku 4 E ) 
» вы даде a: та ws а "НЫ rh у "St 2 4 В ve зи аи "ГЕ UE р S ui ide у р а қ 
та оре Да + EU $ Ы К y 11s E Ы ДЫР Я М қыны Ц ' Te р D NIE | D i a р Қ УЫ | Faris 
са а= РА пие Li RO ; a Ж O: = И ща у LII T a Кола aN E да = A 
— S е ha а м иН Bos Зе аи = 4) "n ЕЛІ TM не. ди ЖИР: | P аа” IP що 
Н b % и Ае | Нана + ae " Е L А Си Ыйы " А П D 
_ ENS _ EA He М я nu s A pu "n. EO "S | 
2. é E L T 4 a, i = 4 А-АА ©) и: `$ “А hs xú Ж y es & < Pax “ғ . ^ г 4 E Й П Й А Ы H . О П А 0 + Й E. 
— Д тре e aum > if я ПА 4 si tese RH k Я e d: О T рғы s m E E mers E ' г i | 
а me RSR аи 5 өрі ue 45, 4: du газа Ты E d ” S Ме Д AE .’ "3 М ory UT A 0 e d ph ы D ui 0 я ; 
ма н ba У A ^ ne a) к % ван "m me = E мр ni ч ком КЕ ge а ы «а J P + Ц О D E m wr Г А 0 P А D Bi 
rt бъда 0 ы ии shé ә? y my Sas 5 жал LE т ^ ~ т? .: 4°: А у Ц НО : О р E y Ц Ры А ы D быр P 
меке É Фа ~. ,* еа ие oon - Boss И е ¢ Е = Й О А А О 
a. мее UE MEOS, pA rd • Ses. E ы Е ЗАЛЕ Ж.Ы ; ." = (7 et e i - 5 М Ы P. й У Ы z у i LE я | ы у и 
OH ме ЕЕ НХ С. пу: aive i E AR it Ка pe D ии. E i А i aunt T E: "v Жы E | | 
Қ "ps См ч 4 T ki „er dá и 5 АЧ К a А На с ^: 20 145% У » BC t х . КЕ ." Й P " X % Ц А П Й 
Vend = a A EQ £ SE t ие DE s "m a ge... T | "EE 
manent сЕ : b = т ей 55 > M 3 M ST E ; E c E D E я к x я Е а р К | at ' В А А, 
н е vi м on m "Mi Хак a А ы “А 1 А er 5 ҰМЫ ~ A heb Я T Se а D а M Й à А , "EU " И А 
— И о ide ay meh > о + eM ; m UN СЕ. E LN TE | MD = , у 
Ка nean oA dL к кла: an ьн да: #75 алалар А Ира y y. Pr 3 И ' ae : PSOE зе A ' А А С 
18 oe ада салса Кери иа E ии ча ME па! es a, ny Т n "DE | ға | u 
4 Sai A едн ~ Ba Е Киа (€ & 19 f8 hs MU da Де, AR ™ <; cx Ы Бао ы T »* s a * НЕ ы Е Ы В: Ц E 
an nn v AN DUE can] a ee 51425 Да ста у що. г. MN "M 
а mnm nd d Ts ыы + Pen eins een ED te e dd Ы q $ Н И l f = у ы}, Ы e . у ы 
а е lod Mh е ur t . қуат өлі 3 *; 2. 4h "d ма M М LM 4. 1 4: р Ни ... E ia: a А Е | у ~ Ы | 
< my T SM Neat: X AE no. ра: 4 и ка Ме Ж | чё. Е па Е. и МЕ и Р я " р » ix у р 22 
E papa аи Li b А г e Да ду ye de MOT pad Pi 4 tum ye es Д аи . Sig s Н ; ' et e н D js * 
с SMS AA MADE зе арта т V E TLN FARM AUR ыы pe И ae . а 0 E " A eae 1 
J carey а ы bi М | РР ыр екн) fn "LI ка ~ Рег јн + UE "T Ц е. “i M E ЫМ "S EH Й , e 7 Е E fi Й а П А В D 
ы, ии Rie ded 4 ARM Pi “4 DM n “~ bon i а қ M Le = en e е Ана .. “e г А Ж MO e Й у М ' и ie d 
ee iy EU: TT пари Ира Be DAE diaeta E АА % ее Га ща M Ae В E О à КЫ > . И и сх ч В: и i” А FE р Ы В A, d 
ENS Ale Sale a ИЕ аа TO E ЫЫ ss e D 456 На па 8: ер М қа E А ' Н ie iat 5 ; 8 . р тож ы у 
ы (54. ar ч ыы nee ey зрее пра Ej ТЕ мр ды ы кызы ы На “ у ne TE e" р у у ie 2 и 0 . 
ЗОНА Se = 2 EE P ; Teg. E LU n А ~ ~ ME E Xe я 
м A тан А , 1 ма “ve aen A 7 ARN ред ü КМ И .. = 5 " "m В E H Й LI M , . d 
e. S T S D m: о Bas E | | d Ж | И ке Е ү 
Ao си а Ее sed m m Ди TUE E "n ши eM | n LN | 
pa SES Ben Deus re de MA ga MA D К le: AS y w 2r у RR КУТ АА . yo” 4 Ы E Áo 7 Es и ы тез 1" > 
ЕЕ пее ры e т е aie iu A E * dk EN грим e M ща > TEE m NC AA 
4 Афера PEEL) a - E Je B L ме м кы, е AN id жм Қы. = Б ЗА А E К M ° Й E Ы " | 
а O ? ELSE |: PON EE UE А. eo пас E ЖЕ: m Зали ще 
a 24 Li и Ы а“. қ un His dis я М %е ANH Е dpe Hee s ша se ы ы ib t е M Pi as .. « .. n 4” “4 ње 4: ' rs d s. EM "E О М М О ы d 2 i 
as th wo ең е у? bates re SA А ME “ ры v ПА ДИ "m E T we [227 ~ . а ~ це M D ГЕС 1 f .. -3 Е Чы n Ы . , О Џ .. • E А у қ 
есет zm о ia. Dos MATES AM ane DUM .` А АМЫ Ыы PT t ae H te ." ДЫ А Й а A , Бы . з ы 1. d Я А " Са a ш ы Е у d А M. | 
ида са Me: nn REIS m RU ORE CER СЫ al o Lo a. DEM dE | 
m H n М % а 1 LI d % B ыы М 7 Ы " че а B = 5 А О Ud В Д 
сас - ОГУ (в mnn AE UE QM TEE p A e n. x | um ME | 
5 abd 1 ~ М | ща б = P m 4 0 та E p P А 6, О Ы 
a] AT EN оран ааа a ка ЕНИ г. че m то пе "ra x | 
ыш e pet 3 soe + Li te * Hee i m " T 17: а m m се: қақ; Hiec de 2 A “ б Е Н ' x Жұм О " М " ПРЕ 2% ы Я А Ы 
BEN о ARA Sun Реза Че Нд di сее: M . 937; ^ d = X Да E , а Уви и 2 E x EL г еи ' у 
— А = EN D eus Un ue TOM. (ера на 220 pro TE. E 
ы + 4 М "n = ^ • О О Ба o . Къ ~ = М е ы ie Ы 
o a br RES n m PR | nm моа i | , | 
ж. Mr — o па E ' у i ' : 
a nn d RUM jen ES b es E ст Em, nmm | | 
Е — k в monte us ER Wu ra doe DE н о Ы m Nu si "an 
ҚЫ Фи а x à з ME A es цени А "ЫЕ 72”. DT .. . E лды d S B ns я 
ze соп ono А Ton OU Hi T d . а t Ki E A P A E : я и" | ка ME I 
- ыы x _ 4; Ж, na "x ie а AX Pu т ұйы Аы ЖЗ OS .” F C. "E у . i ) у | | | 
iy 2 › Ч а 3 IM Jh i 2 M pete] уы 187 E d М P ony г " , Ц . 
мона ша нн Sees E CC tare m Sr Ree ores к y Vat Й ALI ы ее an ie ч s Se ы ыры Ыы р: б "d О Т 4 ili : 2 EE: x ы M : 
ee сее TE UH ES m ADM i: iden Е. пази КА Ери e CE PRG Т” re ү у | | 
rods oos са Baer в Дурин м N (С; posu Pu" T HI о ELE m | m 
cr : on pathy Ero E Е fee ad E SE и m БОК IT um ! 
a пе Ы ае г .0 ЧЕ қын? үөр es р D " пар PP d er ... Др so : Eo ЕН NS ES De 78 ; M n NES. а : 
uL m LI bl (а. Ц ы ч Y r ы E * Й 
c ACE э OE ED Цеци rr INC BE Pu E ZEN КИ : E 
de МА nun и "x e AA ы wand NS RR NM Ива C AM UL DOES AP. PN о En T s пе i P M T | Е у | i a 
x v E Т. AN s HR Toa As у tt TER 52 Ви ЈЕ Ша и н = " JN TM E Er И un Я E ЩЕ А Я ы тағы; d ' 
H P Э rii TM iie 9” Lu ae А па Ки x АМ! DH u N ч Ы v Ыкын nu Ar Ш . АЛ Pur TH ir NM M nds Fr: 215 Je ~ П E Е у и 
AC PAD c оди Бес E zh сака oe e oe А аи a "s DO "nr Dn A x E n E DE ES i 
o. E xen ко за E RA M NE i n Boo ком Га dU ur n HD B | 
| EUN TS ~ QE ^ кор РЕАК ae x S E Дек: ” er м. Я ус АУ E за nv ‚ "n NM E ape що иа р Й Қ | Е : А 
rum 2. DEA po PO EAM но pun T2 | "NGC О в 2 | 
7% а 4 P E АЙЯ ДЕЗ 4 pres 24412 vn рана Wo сащ МЕН e ш PRA Й . \ тыман Ж с Е 4.4 е 
pu nn ыз ie е ERA Kar i ru ие x UN E ааа 2n 
` da mn 15 Г н TM 4 ы A JL) ~ у М M = ` LL v e r ~ Й E . РА Е 
са ma Duae a PS D 2 MEN A б! EARS 2 b Н Я ` ae 7 къ as ae | те | у 4 у 
4 H а Қ 87 à МАЈ om ni А. м ee sw ' ry LT: П “ ise RA ща ur t t M t U [E Г М Ad 1 я ч P2 ы . Ы te Ц 
“м dde d cede: » ecu т я un en eat И Ain y де en ive RA M Meetali Қам си us ve on iu Т р ~ А E ' у : 
ар? f ЖТ LE et Дер MO ene ААН a palit Ка ДД EN B Фа! и рама Mid “1 з © А би w a te тат . A. se E р О у 
^ 3 d + Да AC е. М AA. , Dr СЫ Е. я ees E р » 
ots ”“4 Н ^ з" yay AIV и Нас SER Ма = af NEN ES ^ \, "t и T 4 ы M Дада p. EN я " у я ' Ys у 
ки He а м ¿blas вит носи ки rA QUI UE & mes I ave LED HS T PME . s y E D 
ВО Се porate 2 а пр оси ет Па 2 MES us P 5-55 Em E | 
a Kum т Аз ~ ы) " вест a 28 пин па са v Дада ақыға Mi i POR я t a 50 Т M . А A Р А | g T | 
E EE а Дине си 2 e М uo UNS | Е 9 2 ME | | н 
2-4 3 2 = 4.23 5 b РА ва 4 J d É D D 2 M X П N я D r А 
Дон се: v IN "T 2 Rs УУ а Щи МО ыы г. Rie x 2 | Же; EN "e A en A я 5. ШЕ à : р d 
ie ee 5. Мы bens ГИ Bat т 4 t ГА) "m беди A ВИ Г Tr T yr EET Rol ui ми и Дам " Ра у "9. con Pe erii а, КШ ` б ¿> EI А Й 
ма MP hind е“ bi РАМА: MA AS Foe по el А T Ae [ Ud Ы; A E TR uie Ut ML e за КН б; ГМ: : Ы M i 
eu + n nr er 2. Md bit 1] СТ НЫ А neu m S" ek дето 5. А "bd. ки ҚЫЗ, Мы B Қы E Б M. А | А А i " б М E "E: В 
Aft RIP В ЕА one Ar 4 "rti Ku Ех е M3 у» Ar ner А Ке А vi ha m асасы PEN Ди 4 Ар js КС as rn E we A n А "e О А Ш 0 
el и чада ДА n и пра i ҚАҚ: YN un T A + a an wae ғ P СЕХ DM Te Er бу DX қ и | ' : 
соса 557% ХЕ ad Ка d eA PUN FIM X НЕ? А aes En d ща d S s o ru Ч ' L ME б. х ; ; 
с por ЖЕ PRA EIN RECTE es Bodo Ar и НА DOLO Ди ў | зар lp AEN : AS н = 
ғ Иа Юри си is "n TEN А. TA Nn к ÑO УР M р n вт ТЫ” " aĵ 0 LN ~ Де t, М АЫ Е EET a b г Й ы ' HAE a Es А А 
ККД ЖЕНУ: Gee Аы ' ране v Да Қ ЈА M ди бор: пей N 5. .. да Ju я. ГИД cu EL i x Ac n m а А D LM ha 
X K 78 ^ | Le Wes м П 2 . . б " А Е, Ы H 
mora Е ТА ДУ е пи ND а P DO E UD DD x a аи nt TER QUU on | 
Lag 9? VES Ру n ЗА Қ E Y ыш A e КА Ы e en и ' f Ааа Ц МАЈ EM па қы теи E З ' ща. E у ra m м: в 
М m А Қ. А v Жж А eE М Bin Би “N DM П жы " uM ҚОЖ, " H П Й LI А d PU қ А N А . 
дос ЊЕ Y 44:44 "m Ы ЙЫ: Dif САН мс ee ме" D KENA r He D ШЫ $. : Б, E D О is .. у 0 ы Ш P 
ШАРА. MES Od ЖЕТЕС ха “+ .7 ' коса Di ael "ue ‘ b " Ш " A. EN ' . 
МА | m # et ALL С et Да 7 cr NS A " AL О В ы О Ы Џ k 
ди | а РА ЕЕ гр 78 E "UP У П а PUR. ДЫРЫ N ' АТИ H d d D В „ , . 7 А А 
D fr Ae Aeg 19" v т! TAE vd ue Џ И : PL a! ез о y О ы И i Ы р А у 
ъф fi d dn 4 T Је 441 ete al + ‘ \ M $t us Ж. АЦ L ye Ы Ц E М ji, ie " D ' d ` PM Е 
НА "n \ вые Miet мы | г T HE * ДР зи ДА to 0 m ‘ E қ А Т и pe . " А Қы; Ш 
К we ^y A Nx, Vn 34 s M. БИ ту o a P 4% hr) | S Т пи ' ы n 
е S 4 PY КЫ is Оу КЫ, М E 2 (C sas à^ U a 3 "T е шы К ы П | de d Џ б E Е Ц . М ME " 
ULA > “к! ни t Erde P , = ы К . • D 0 $ А D О О а О 
a XT 54. | А EO ID Бру COD А er M M P PA А [ET] e Аш ws Па | ПЕ A ч А à L | 
ы Mae n m 2 Ну HI e ó A Р) be 5 , A, ( ү M қ Y ж М ei. T за f М " r р ти Й ida Й Ы 
H y Far E 9 2 ww күр ы ier М m Ц Ы , 1 ы 1 .. А A Ц d У ы 0 
Ps С "m Nur « 4 ' ok й D D "n EP i ' Й ru 
MT ! На Мен М oS 7 uuo 321 à 4 e я Ы 0 ' ща, 
S Y DE s Му > 1 ти 4 MU Й Фен "n AM a 41 З N Ds E Е А П b. 4 SU " П О 
rt А na ВОИНА Ted 945 Е TAM Se ж” ^ i р i р Д у р : р d 
КЕ er it re 4 м.п А > ы A я DL D ы е та a D ғы se О у 
ЖИ И £3 DUE rim Ls " ih st T ИГЕ | D D n 
ЕТ sn EUM UAE "T Ae is " А 5 Ц .. ” Г А ШЫ Ы ' 
Uu EN Џ | Fr n О Ги 24 "ie "n . G А mu И А ^ n E Ш А = И 
.4 ae АП n ба Т " . . а ri і ' М E s Й Ы "Li К | 
ery te a " A M P 10 и " П Я ЕРИ А Р " ' 
е 4.» t uM П EL ы ғы А0 ra П ' р v 5 P ш y mm А В , 
Қ $ г г. r ie LI П Hans: и ' А y Џ d 
» H 1 = ПАДА, А , РГ Ед: 4 А yf Р 
TAE Е ТА m Е ? Й f р , TEE LIN , A) 
ve PII HR ща за А ара П НЕ " А П net E 0 + Е Ы А EJ 
LI ba GT : Sur (t и us АІ; гуа. P я р А | 
А || „у КА ЗА ДА Н АА x "E Ны > 4, П О 
Е ды A Qd ' ' DR Й D А 
З AL E Pl Да E К я P ES У Џ " D z 
N n Ы "m E "m. А " 
е з М ur: Ва s Far A 2 ' Ж E : AR L ' | | 
. Er А Я р А ELE П А ; | | ' О 
е ЈА ta Peg LE у d Й у r Ц 
AM Js e "t П 4 Ы = Ду О | Е | 
М. . LI H- 
ы AS ~ ^ р 4 ' "n О 5 қ 
LP d dx М ba il. 1 . П 
. 4 Aa К ҚЫ з: Ри у Ы 
ДЫ Ва ы А D 
LI Ц t 
LI 1. Ll + 
. ы ы 5 n t: П PLI 
Н ЕА, ЧЕ р 
"E : 
( n! О , р 4 М 
D ss О ы ' 
р | | 
ҒЫ М Ru y > 
b П О 
ка М * 
р 
G Џ 
Џ | n g 
L] 
' 
L] 
'" у 
! 





